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Translation Editor’s Foreword 


The techniques of many-body theory are best known in their 
application to the electronic states of metals and to statistical 
problems. In the theory of heavily doped semiconductors both 
the quantum and statistical mechanical aspects of the problem 
are prominent. The present tract is a review by Professor 
Bonch-Bruyevich, one of Russia’s leading solid state theorists, 
of the quantitative progress which has been made in this chal- 
lenging field. The English edition is a translation of the Soviet 
original which appeared in 1965 and should be of value as a 
reference to all physicists working in this field. 


R,S. KNOX 
ROCHESTER, N.Y., MAY I966 


Preface 


The current interest in solid state theory is, in the main, 
the result of confluence of two factors: recent progress in the 
problem of many bodies in quantum-statistical mechanics, and 
demands of practical semiconductor physics occasioned, in 
turn, by the ever increasing use of solid state devices in sci- 
ence and technology. 

Professor Bonch-Bruyevich has rendered a real and well- 
timed service in pulling together the available information on 
heavily doped semiconductors, because the problems associated 
with these materials are among the most pressing in modern 
applied and theoretical physics. The reader will find here a 
detailed review of the many-electron approach to the theory of 
doped semiconductors, as well as a systematic exposition of 
the many studies done in this area (much of it by the author and 
his students). Some of the interesting topics covered are the 
energy spectrum of heavily doped semiconductors, pulse dis- 
tributions in them, as well as screening effects and absorption 
of electromagnetic waves. The necessary mathematical ma- 
chinery of the theory is explained simply and concisely. 

It is hoped that the book will prove of value not only to 
specialists in the field, but also to a wider group of solid-state 
theoreticians and experimentalists. 


S. TYABLIKOV 
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CHAPTER 1 


Introduction 


The study of the electronic properties of semiconductors 
containing a rather large amount of impurity has become 
extremely important in recent years.* This is primarily as- 
sociated with the discovery of tunnel diodes (see, for example, 
the review [1]); however, thearea ofapplication of heavily doped 
materials is actually much broader. These semiconductors are 
used in, among other things, lasers and thermoelectric devices, 
semiconductor catalysts and, of course, thin films. The term 
‘‘impurity’’ is not necessarily to be interpreted literally—it can 
refer to any structural point defects that lead to a disturbance 
of strict periodicity of the lattice. Disordered alloys (including 
metallic materials) represent the limiting case of materials of 
the type under consideration.** 

Obviously the study of heavily doped semiconductors is a 
part of the more general problem of disordered structures. The 
study of the latter is of nearly primary interest in present-day 
solid state physics. 

The phonon [2-5] and electron {6] spectra of materials of 
this type have been studied for the case of relatively low im- 
purity concentration. The method employed in the works cited 
represents, in essence, a reasonable perturbation theory based 


*The exact meaning of the expressions “rather a large amount of impurity,” “heavily 
doped semiconductor,” etc., will be explained later. 

*"In essence, any impure specimen (the only thing semiconductor physics dealt with 
in its early stages of development) is heavily doped. However, serious research was 
begun only after it became possible to obtain materials doped with a specific impurity 
and at a specific concentration. 
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on the expansion of certain quantities in a series in powers of 
the impurity concentration. One or another of the characteristics 
of an ideal crystal, assumed to be known, is taken as the zero- 
order approximation. In this fashion it has been possible to ex- 
plain a number of the important properties of slightly disordered 
systems and to obtain a quantitative description of their energy 
eigenvalue Spectra. However, in the caseoftheelectron spectra 
the conditions for the applicability of the series expansions are 
somewhat more restrictive than one would like for consideration 
of most materials of experimental interest. Thus itwill be seen 
later that in the semiconductors employed in tunnel diodes and 
lasers (Ge, GaAs, InSb) the impurity concentration is greater 
than that required for use of the expansions, In this instance the 
doping element is more of a microimpurity so that a statement 
of the problem based on the specific properties of the correspond- 
ing ideal crystal makes complete sense. Moreover, some specific 
effects that require special study alsoappear. Theseeffects will 
be examined in the next three sections; a definition of the con- 
cept of a “‘heavily doped semiconductor’? will also be given. 

Later (Chaps. 5, 6) a mathematical device suitable for de- 
scribing the systems of interest to us will be considered. Within 
the limits of the present review only a discussion of the type 
mentioned above will be given; a more detailed knowledge of the 
device can be obtained from the papers cited in the appropriate 
chapters. The reader who is interested only in the results of 
the theory can omit Chap. 6,aswellas all the subSequent mathe- 
matical portions, limiting himself only to the summaries of the 
results located at the end of Chaps, 7, 9-15. Chapters 7-11 and 
Chap. 14 are devoted toa study of the energy spectrum of heavily 
doped semiconductors and the associated problems. Some prob- 
lems in kinetics which are important for an interpretation of 
many of the experimental data are considered in Chaps, 12, 13. 
Finally, in Chap. 15 the characteristics of intrinsically inhomo- 
geneous systems—the contact region of two heavily doped semi- 
conductors with identical or different types of conductivity—are 
investigated. 

It should be noted that the theory of heavily doped semicon- 
ductors in its complete form started to evolve quite recently 
and many important problems have yet to be solved. Thus, for 
example, galvanomagnetic effects have not been completely in- 
vestigated from a theoretical point of view. Explicit equations 
have recently been derived for the Hall and thermoelectric 
power constants of a heavily doped semiconductor [118]. The 
corrections to the equations of the standard theory are small 
in the limiting case of heavy doping. 
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Some of the usual concepts of semiconductor physics, used 
consistently for interpreting experimental data, are inadequate 
for the problem in heavily doped materials and sometimes simply 
become meaningless. The need to call this circumstance to the 
attention of researchers, as well as the desire to draw the at- 
tention of theorists to the complex and interesting problems of 
disordered systems, served as the primary reason for compil-~ 
ing the present review. 


CHAPTER 2 


Impurity Levels and 


Impurity Bands 


It is simpler to approach the definition of a ‘‘heavily doped 
Semiconductor’’ by considering the change in the band spectrum 
of the crystal as the impurity concentration increases. It is 
well known that isolated impurity atoms can lead to the creation 
of discrete levels located in the forbidden band. The wave func- 
tion of an electron occupying such a level is localized in the 
vicinity of the impurity atom. The energy spectrum of the 
semiconductor in this case of ‘“‘light doping’? is represented 
schematically in Fig. 1 (the ‘‘allowed’’ regions ofthe continuous 
spectrum are crosshatched). For definiteness let us assume 
that each impurity atom gives rise to only one discrete level. 
The symbols are the usual ones: E, —the bottom of the conduction 


wt 


Fig. 1. Lightly doped semi- 
conductor. 
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band, E,»—the top of the valence (i.e., the bottom of the forbidden) 
band, E:—impurity level. Obviously the concept of impurity 
atoms isolated from each other can be valid only as long as the 
‘‘repions of localization’’ of the associated electrons do not 
overlap to an appreciable extent. (For this to be true the im- 
purity concentration n must be quite small.) When such over- 
lapping sets in, the impurity levels should ‘‘spread’’ in the same 
way as do the levels of the atoms of a crystal relative to their 
discrete nature in the vapor phase.* 

A complete band of allowed states (separated, however, from 
the conduction band by a forbidden region) appears in the for- 
bidden band. It is commonly called the impurity band although 
it must be admitted that this is an unfortunate designation: 
nothing is said here concerning the dependence of the energy on 
crystal momentum. Moreover, using the concept of crystal 
momentum here requires justification because, generally speak- 
ing, the impurity atoms are not periodically arranged. Con- 
versely, it is more reasonable toassume that their arrangement 
in the lattice is random. In a macroscopically homogeneous 
crystal the average impurity concentration is, of course, con- 
stant over a relatively large region; however, insmall volumes 
(for example, in a sphere with aradius of the order of the radius 
of the effective interatomic force) the variation in concentration 
can be considerable. This circumstance leads to a second 
mechanism for the ‘‘spreading’’? of levels. That is, we isolate 
some arbitrary impurity atom and consider the effect of all the 
other impurity atoms on the electron localized near it. A vari- 
ation of concentration as we consider different arbitrary atoms 
obviously leads to a variation in the potential energy of the as- 
sociated electrons, i.e., ultimately to a spread in the ionization 
energy. This may be detected, sometimes optically, as a 
broadening of thelevel. Todifferentiate it from the first (‘‘quasi- 
chemical’’) mechanism for the development of an impurity band, 
this second mechanism can be called the ‘‘van der Waals”? 
mechanism, 

We will designate as moderately doped the semiconductor in 
which there is an impurity band separated by a forbidden region 
from an “‘intrinsic’’? band in the crystal (the conduction or 
valence band). The energy spectrum for this situation is shown 
in Fig, 2 (the region corresponding to the impurity band is 


*This analogy is not complete because the fact that the impurity atoms are located 
aperiodically is ignored. Moreover, the assumption of a periodic distribution of the 
impurity atoms (with a period equal to a multiple of the lattice periodicity) would im- 
ply the absence of impurity scattering, which is clearly at variance with experiment. 
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Fig. 2. Moderately doped 
semiconductor. 


crosshatched). Drawing sharp edges on the impurity band is 
somewhat conditional, they, too, can be smeared out. 

The impurity concentration at which the development of an 
impurity band becomes noticeable depends on the nature of the 
atoms and can be calculated only within the framework of some 
specific model. One of the methods for experimentally detecting 
impurity bands involves a study of the low-temperature electrical 
properties of the crystal. When an impurity band exists, the 
electrical conductivity does not become zero upon ‘“‘sweeping 
out?’ electrons from the conduction band (or holes from the 
valence band) but remains finite, though small (see, for exam- 
ple, {7]). The Hall constant R also exhibits a distinctive tem- 
perature dependence. As is easily verified, in the presence 
of two bands with carriers of only one type it is given by the 
equation [7] 


_ AyWyR yy thee ory 24 
a ec (yy +2 p12)? ( ° 


Here m, M1, Mi, and no, pe, fe, are the concentration, drift and 
Hall mobility of carriers in bands 1 and 2. If the temperature 
dependence of the mobility is neglected, Eq. (2.1) has a maxi- 
mum when 


Nyy =Nepe. (2.2) 


At the same time, in the absence of an impurity band (y2.=0) we 
would have had R~'!/n,, i.e., the Hall constant would increase 
monotonically with a decrease in temperature. 
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Using the above reasoning it was established that in the 
particular case of germanium doped with group V donors a 
spreading of the level becomes noticeable at an impurity con- 
centration of n ~ 10'5 cm-3,* 

The ‘‘width’? of the impurity band obviously increases as n 
is increased. This can be observed experimentally as, for ex- 
ample, a decrease in the activation energy which determines 
the concentration of electrons in the conduction band [8]. Ata 
sufficiently high impurity concentration the activation energy 
finally becomes zero. (In germanium and silicon, doped with 
group V donors, this occurs at n ~ 3x10! cm-'andn ~ 10'8 
cm’, respectively {[8, 9]). This would indicate that the upper 
edge of the impurity band reaches into the conduction band 
(Fig. 3). Under these conditions, strictly speaking, one cannot 
speak of an impurity band and of a conduction band literally—a 
single continuous spectrum emerges. We will designate as 
heavily doped the semiconductor in which the impurity band 
overlaps the intrinsic band next to it. 


Eo 


Fig. 3. Heavily doped n-type 
semiconductor. 


The first important characteristic of heavily (and moderately) 
doped semiconductors follows directly from the definition: the 
impurity simultaneously shapes the characteristics of the charge 
carrier spectrum and causes their scattering. We note that in 
the opposite case of light doping the electrons localized at the 


*Here we are notconcerned with the question of the conduction mechanism in terms 
of the impurity band. lt can vary, depending, in particular, on the mechansism by which 
the band is formed. For our purposes it is sufficient to bear in mind that any factor 
that provides a possibility for direct electron transitions between impurity atoms also 
leads simultaneously to a broadening of the level. 
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impurity levels do not take part in the migration process, and 
the structure of the conduction band is practically independent 
of the presence of impurities. In other words, in the case of 
low impurity concentration the impurity atoms shape only the 
spectrum of the bound charge carriers and cause scattering of 
the free carriers with practically no change in the latters’ 
energy spectrum. In the case of heavy and moderate doping the 
impurity plays both roles simultaneously. 

A second characteristic of heavily doped semiconductors is 
also evident: in the absence of compensation a degeneracy of 
the electron gas, even though partial, should be expected.* 


*Here and in analogous cases later on, we will consider the material to be n-type 
for definiteness. All considerations carry over to the case of holes by a simple change 
of words and symbols. 


CHAPTER 3 


Neutralization in Terms of Impurity 


Configurations. Crystal Momentum 


A third characteristic of heavily (and moderately) doped 
semiconductors can be exposed by noting that most of the elec- 
tronic characteristics of the system that are directly calculated 
from quantum mechanics depend on the arrangement of the im- 
purity atoms. The reason for this is obvious: in the solution of 
the appropriate dynamic problem the coordinates of the impurity 
atoms can mostly be assumed to be formally given because they 
also enter the answer as parameters.* On the other hand, the 
actual coordinates of the impurity atoms are certainly not given 
in an experiment and the ‘‘answer,’’? containing Ssomany unknown 
parameters (~ 10!8§, for example) is of no interest. It should 
also be remembered that in different parts of the crystal the 
impurity can be distributed differently, whereas the result of 
a macroscopic experiment refers to the entire crystal and not 
to individual small portions of it. In other words, an averaging 
of the desired quantities over the entire volume of the crystal 
occurs in the macroscopic experiment. The third important 
characteristic of heavily (and moderately) doped semiconductors 
consists of the requirement of this averaging. 

Let the total number of impurity atoms in the crystal be N, 
and their radius vectors be R,,...,Ry. Let us designate by 


*In the case of light doping one can solve the dynamic problem concerning the be- 
havior of electrons in the presence of only one impurity atom and then “multiply the 
resulting value” by n. For this reason the question of the location of the impurity 
does not appear—one can simply assume that the given impurity atom is located at 
the origin of the coordinates. 
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P(R,.... Ry) dRy...dRy the probability that atom 1 is located in the 
volume element dR, about the point R,, atom 2 in the volume 
element dR, about the point R2, etc. Then the average value 
(over the impurity configurations) of any physical quantity A(R, 
.., Ry is defined by 


A= (4. ...,Ry) P(Ri,...,Ry)dRi...dRy, (3.1) 


where the integral is over the entire volume 2 of the crystal. 
Obviously the normalization condition 


\o PRL... ,Ry)dRi...dRy=! (3.2) 


must be imposed. 

In what follows (with the exception of Chap. 15) we will as- 
sume that the impurity atoms are arranged uniformly on the 
average, in the lattice. This means that integrals of the type in 
Eq. (3.1) are not taken over the entire volume of the crystal but 
only over some relatively large portion of it that is independent 
of the location of this portion in the lattice. A lattice for which 
this method is valid is said to be macroscopically homogeneous. 
Its symmetry is obviously the symmetry of the corresponding 
ideal lattice. This circumstance allows one to introduce the 
concept of the crystal momentum of the charge carrier. Of 
course, because of the translational symmetry of the lattice 
the one-particle density matrix, p(r, r’), has the following 
properties: 


p(r, r’) = p(r +a, r’) = p(t, r’ +a), (3,3) 


where a is a lattice vector. Consequently, the eigenfunctions of 
o(r, r’) are Bloch functions, characterized in particular by the 
crystal momentum p (see [10]). The equation for these eigen- 
functions in the coordinate representation has the form 


(lr, ) ple’) de’ = Ep (r, (3.4) 


where € are the eigenvalues. Making the substitution r+r+a and 
bearing Eq. (3.3) in mind, we see thatthe function (r-+a) satis- 
fies the same equation with the same boundary conditions. The 
eigenvalues &, are obviously characterized by the crystal mo- 
mentum p and the ‘‘band index’? !; as is known, in a system of 
Fermi particles they are equal to zero or unity. It is simple to 
verify this by employing the well-known representation of p(r, r’) 
in terms of the particle creation and annihilation operators @ 
anda ({11, 12]}): 
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p(r, r’)= (a(ra(r’)>. (3.5) 


Here the symbol <...> designates an averaging with the total 
density operator: 


¢...) =Spur(p...), (3.6) 


where the spur is taken with respect to second quantized states. 
In order to obtain the intrinsic ensemble, let us consider the 
case when only quantum mechanical averaging remains in Eq, 
(3.6) [as occurs, for example, under conditions of thermo- 
dynamic equilibrium when the temperature of the system ap- 
proaches zero]. Changing over to the representation in which 
the matrix p(r, r’) is diagonal, we obtain from Eq, (8.5) 


c= Ca 1Ap I >om=Mp I- (3.5) 


Here the index ‘‘QM’’ on the average symbol means that one 
must understand the average only in the quantum mechanical 
(but not in the thermodynamic) sense; fp is the number of 
particles in the state(p, /). In a system of Fermi particles this 
is zero or unity. 

The system of eigenfunctions of Eq. (3.4) is especially con- 
venient for computing the average values of physical quantities 
represented by additive operators (for example, the current 
density). Actually, in accordance with the general rule for 
computing averages [11], we obtain for the average current 
density 


<j>= 2 DO, Nile 1) < apap >p, (3.7) 
p,p’ 2,0’ 

where (p, /|j|p’, /’) are the matrixelements of the current-density 
operator in a system of Bloch functions, and the index ‘‘f’’ on 
the average symbol means that one must perform the averaging 
operation while taking account of external electric and magnetic 
fields. The single particle density, <a,-)-a,.>;, generally speak- 
ing, is not diagonal in the (p, )) representation. In weak electric 
fields, however, it can be expanded in a series in terms of the 
field intensity [12]; then Eq. (8.7) (in the absence of a magnetic 
field) ultimately assumes the obvious form 


>= Di, (lilp, Dar. 9, (3,7') 
pl 


where n;(p, !) is the average number of particles in the state 
(p, l), calculated while taking account of the external fields. 
Equation (3.7') justifies the name “‘charge carriers’? (or free 
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charge carriers) frequently given to particles in states of given 
crystal moments. 

An equation like (3.7) is also valid for any other additive 
physical quantity. Let its corresponding operator be A. In the 
absence of external fields the exact formula is clearly 


(A> = 3p, L[A|p, 4) <aprapr>. (3.7°) 


pil 


Here, the sum over p, as in Eqs. (3.7) and (3.7), is taken over 
all points of the Brillouin zone allowed by (e.g.,) periodic bound- 
ary conditions at the faces of the crystal. 

In an ideal crystal the Bloch functions , satisfying Eq. 
(3.4), are simultaneously eigenfunctions of a Hamiltonian (in 
the one-electron or ‘‘quasi-particle’’ approximation). Conse- 
quently, one can assign a specific energy to the charge carrier 
and introduce the dispersion relation in terms of the dependence 
of the energy on crystal momentum. In a semiconductor with an 
impurity, however, this is not so: the Hamiltonian does not 
have translational invariance (the averaging procedure has 
meaning only in connection with observed quantities, not the 
operators). In other words, unlike the case of an ideal crystal, 
states with a given crystal momentum in a semiconductor with 
an impurity are not strictly stationary. Scattering at impurities 
leads to their attenuation. Consequently, the dispersion law in 
terms of p has only an approximate meaning. In itself this cir- 
cumstance is quite trivial, since, in an actual crystal, states 
with a given crystal momentum are not always stationary, but 
at most quasi~stationary. However, it becomes significant and 
important if the first characteristic of the systems under con- 
sideration mentioned above is borne in mind: it is impossible 
to consider the band structure of a heavily doped semiconductor 
without taking account of the interaction of the charge carriers 
with the impurity. However, introducing this interaction leads 
not only to a change in the energy spectrum but also to attenua- 
tion. As we will see in Chap. 7, the attenuation constant may be 
by no means small compared with the corrections to the spec- 
trum. In this sense a heavily doped semiconductor could be 
called a medium without a dispersion relation.* This is its 
fourth important characteristic. 


*As we see in Chap.8, in the limiting case of very heavy doping the attenuation in 
the energy region close to the Fermi energy is, however, quite small for many prob- 
lems in kinetics. In this situation one can speak of a spectrum but it is almost the 
Same as in an ideal crystal—the corrections due to impurities are also small in this 
region. 


CHAPTER 4 


SCREENING 


An increase in the impurity concentration is accompanied, 
as arule, by an increase in the number ~ of free charge carri- 
ers. Consequently, the interaction between them can turn out to 
be significant. One of the effects of the interaction—screening 
of the field of the impurity atoms by free charges—plays an 
especially important role. Its nature is simple to understand by 
considering a somewhat special (although experimentally feasi- 
ble) system, namely, a lightly doped semiconductor under high 
injection conditions. In other words, let us consider the prob- 
lem of the energy spectrum of a homopolar semiconductor con- 
taining one donor and many electrons in the conduction band; 
system neutrality is ensured either by a corresponding number 
of holes (injection), or by other donors which, in accordance 
with the condition of light doping, cannot be explicitly taken into 
account, 

We are restricted to the case of ‘‘shallow’’? donors (such as 
group V elements in germanium or silicon). Then the effective 
mass method [28-32] can be used, and, for the sake of simplicity, 
we will assume that the effective mass tensor is isotropic (this 
limitation is not fundamental for our problem). Then the energy 
spectrum of the charge carriers is determined from the 
Schrodinger equation, 


— = VptV () p=Ey. (4.1) 


Here m is the effective mass, V(r) is the potential energy of the 
interaction of an electron with a donor located at the origin of 
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coordinates, and the EF the energy eigenvalues measured from 
the bottom of the conduction band. The region E> 0 is the 
conduction band; negative eigenvalues, if they exist, form a 
discrete set and describe the impurity levels in the forbidden 
band. The function V(r) is the potential energy of the interaction 
of a carrier with an ionized donor screened by other free charge 
carriers.* In the absence of the latter (and, therefore, in the 
absence of screening) it could be approximated by the Coulomb 
equation, — e*/er, where « is the static dielectric constant of 
the lattice. When account is taken of the screening effect, the 
potential of the donor becomes short-range, Its explicit form 
has been calculated repeatedly by various authors in connection 
with different experimental conditions. In particular, in a non- 
degenerate system the general Debye equation [13] gives a 
satisfactory approximation: 


Vin= — exp a) (4,2) 


lo 


where ro is the Debye radius, ry)= (eT/4xn,e2)'*, and 7 is the 
temperature in energy units. 
For our purposes it is sufficient to know only that when 


\arV(r)|< oo, AS 


is known (see, for example, [14]), under such circumstances 
the number of discrete eigenvalues of Eq. (4.1) is always finite 
(in contrast to the purely Coulombic problem); in fact, for a 
sufficiently rapid decrease in |V(r); as one goes out from the 
origin this number is equal to zero. The rate of decrease of | V| 
is determined, obviously, by the concentration of screening 
carriers and the temperature; in particular, in anondegenerate 
system r,>~n,'” in (4.2). The physical meaning of what has been 
said is that screening reduces the interaction ofan electron with 
the donor, as a result of which thelevels of the latter, when the 
concentration of charge carriers increases, are ‘‘pulled closer”? 
to the conduction band and, finally, disappear [15, 16]. From 
physical considerations it is evident that the disappearance of 
the levels takes place when the screening radius ry (or some 
other characteristic length which determines the rate of de- 
crease of |V| in the case ofnon-Debye screening) becomes equal 


r—-oo, |V(r)| decreases faster than r™’, i.e., 


*Strictly speaking, it 1s not clear a priori to what degree it is valid to take account 
of the essentially “many-electron” screening effect in the formal ‘‘one electron” equa- 
tion (3.1). This problem will be discussed inChap. 5. For the present, however, we 
must depend on intuitive considerations in favor of such a possibility. 
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to the first Bohr orbit in thecrystal, ay) = eh?/me?. For the Debye 
relationship of Eq. (4.2) an exact calculation was carried out on 
an electronic computer [17]. In a nondegenerate system for ¢ = 
16 and m = 0.2 mp (mo is the true electron mass) the kinetic con- 
centration of carriers at which all the discrete levels disappear 
is about 10" cm’. Approximately the same numbers are obtained 
in the case of a degenerate gas of screening charge carriers. 
It is clear from the previous section that this is just the con- 
centration that characterizes a heavily doped semiconductor (in 
every case, in the absence of appreciable compensation). There- 
fore the question of the role of screening in heavily doped 
semiconductors acquires foremost significance. 

It is evident that under the prevailing conditions for level 
disappearance it is pointless to discuss the problem of impurity 
bands because of their broadening. From this, however, it does 
not follow that the concept of impurity bands is not general. In 
the case of a heavily doped semiconductor the problem differs 
only in that the charge carriers should not be supplied externally 
but by the impurity itself. If it is assumed beforehand that all 
the impurity atoms are ionized, then the screening actually leads 
to the disappearance of the donor levels. If, however, it is as- 
sumed that all the electrons are localized, each to its ‘‘own’’ 
donor, then there will in general be no free charges—and, 
hence, no screening. This means [16] that in any given case the 
problem demands a self-consistent approach: the solution of the 
question of the existence or absence of impurity states depends 
on the distribution of charge carriers, which, in turn, is itself 
determined by the form of the energy spectrum of the semicon- 
ductor. In other words, in solving the problem of the energy 
spectrum of a heavily doped semiconductor it is impossible to 
make a priori assumptions concerning the spectral character— 
whether impurity states exist should be derived from a solution 
of the dynamical problem of the system ofelectrons and donors. 
It is evident that discrete levels are not involved in any self- 
consistent solution of the problem in the case of heavy doping; 
the overlapping of the corresponding wave functions when n~ 
10!7 - 10'§ cm-* would not be small, i.e., a state with a highly 
localized electron at some impurity atom is neither stationary 
nor quasi-sStationary. 

The requirement of a self-consistent (in the sense stated 
above) and, consequently, fundamentally many-electron state- 
ment of the problem comprises the fifth important character- 
istic of a heavily doped semiconductor. 


CHAPTER 5 


Characteristics of a Many- 
Electron System 


AS we saw in the previous section, the study of the energy 
spectrum of a heavily doped semiconductor is a problem in the 
many-electron theory of the solid state. Let us consider here 
some of the basic concepts required for what follows. Only a 
Summary of the results used in the following sections is set 
forth, Proofs of the statements which will be made here can be 
found in the literature cited below. In particular, the words ‘“‘it 
can be shown’’ without a source citation refer to [12]. 

First of all it is necessary to define correctly the concept 
of ‘‘charge carrier.’’ Within the framework of the one-electron 
approximation this is trivial—the charge carriers are simply 
the electrons in the unfilled bands. However, an interpretation 
is required in the many-electron statement of the problem. 
This has been given by various authors [18-22, 12] from points 
of view that are formally somewhat different butare substantially 
equivalent. 

The concept of the solution is suggested by Eqs. (3.7)-(3.7 ), 
and also by the circumstance that in an intrinsic semiconductor 
at low temperatures the charge carriers disappear. It is nec- 
essary to consider the charge carriers as elementary excita- 
tions which are characterized by a specific crystal momentum 
and, possibly, by some kind of discrete quantum numbers. In 
other words, one must ascertain how the energy of the entire 
many-electron system is altered through a variation of one of 
the numbers n(p,/) about unity. It is natural to designate this 
change in energy, W(p, J), the charge carrier energy; how- 
ever, this name requires justification: by ‘‘energy’’ one usually 
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understands a characteristic of a stationary state. Consequently, 
it is necessary, after having assumed that the State (p, /) is 
produced at the instant of time ‘=0, to see how the probability 
of detecting it at time ¢>0 will vary with time. The amplitude of 
this probability is clearly 


Apu (f) =< Qpz (£) ap: (0). (5.1) 


From general considerations it is evident that this quantity will 
either oscillate while remaining constant in absolute value, or 
decrease with time, In the first case the State(p, /) is stationary 
and the concept of charge-carrier energy has a precise mean- 
ing.* The statistical properties of the system of carriers are 
no different from the properties of an ideal Fermi gas, In the 
second case the State (p, /) is not stationary and, strictly speak- 
ing, it is impossible to introduce the concept of charge-carrier 
energy. It can have only an approximate meaning, if the attenu- 
ation is sufficiently small, and the ‘‘level width’’ must be small 
compared to the energy W(p, /) measured from its minimum 
energy (from the ‘‘bottom of the band’’). Henceforth, when 
speaking of the charge-carrier energy spectrum, we will be 
referring to its attenuation constant as well as its energy. 

In a semiconductor the procedure can be divided into two 
Stages. In the first stage in an intrinsic semiconductor the 
elementary excitations of the many-electron system—conduction 
electrons (called simply electrons from here on) and holes—are 
distinguished.** The minimum energy, A, necessary for form- 
ing an ‘‘electron-hole pair’’ is that which iscalled the forbidden 
band gap in the normal band diagram. Strictly speaking, the 
number of electrons and holes is not a constant of the motion. 
Pairs can be created and can disappear as a result of fluctua- 
tions in the force field in the crystal through thermal motion 
(we now retreat from every type of external effect). However, 
for sufficiently large values of A (typical for semiconductors) 
the probabilities of such processes are extremely small.*** 
Consequently, in many problems the number of electrons and 
holes can be considered, in practice, as a constant of the mo- 
tion. Thus the problem is reduced to the study of a mixture of 
two gases, namely, conduction electrons and holes. These gases 


*lt is still necessary to prove that the oscillation frequency is W/A. This proof 
is easy to do, for example, by the Green’s function method. 

**There are still exciton-type excitations, but they are of no interest to us here. 
***This fact is immediately appreciated from experiment: it is well known that re- 
combination times, as a rule, are extremely small compared with #/A(A™ 1 € V). 
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are nonideal because the particles comprising these gases in- 
teract with each other in accordance with Coulomb’s law. The 
effect of the lattice is manifested in a) the periodic force field, 
acting on the electrons and holes, and b) the attenuation of the 
Coulomb interaction between them because of polarization of 
the medium, 

In a homopolar crystal effect b) can be described approxi- 
mately by introducing the static dielectric constant e of the 
lattice into Coulomb’s law. This approximation is validuntil the 
frequency W/h, where W is the energy of the electrons and 
holes, each measured from the minimum of the corresponding 
zone, lies sufficiently close to the dispersion region of «. In 
other words, the excitation energies must be sufficiently small. 
For most electrons and holes in germanium-type semiconductors 
this condition is usually satisfied. 

In a nonideal semiconductor the segregation of the elementary 
excitations, the total number of which is practically unchanged, 
actually proceeds in this manner. A trivial difference arises 
only in the fact that now, for example, pairs of ‘‘conduction 
electrons and holes at a donor’’ can play the role of elementary 
excitations. Let us return to the problem of a nonideal gas 
whose particles interact both with each other and also with the 
impurity atoms (the latter interaction also being screened by 
the static dielectric constant ¢). 

The study of the nonideal electron gas (or holes or both) 
constitutes the second stage in the solution ofthe problem of the 
Semiconductor energy spectrum. For definiteness let us con- 
sider an n-type material, disregarding the presence of free 
holes (a one-component gas). As can be shown, the charge 
carrier energy spectrum is determined by the eigenvalues £ of 
the ‘‘effective Schrodinger equation,”’ 


{E+ yt— UL) -U() bo — 


= \dr'M (r, r’, E)p(r’) = 0. (5.2) 


Here Ux, and U are the potential energies of an electron in the 
periodic field of the lattice and in the screened field of static 
defects (for example, impurity atoms); M (r, r’, £E) are the matrix 
elements of the mass operator describing the interactions of 
the electrons with each other and also, possibly, with phonons; 
y iS an auxiliary function (under certain conditions [12] it can 
be used as the true wave function in calculating the average values 
of additive physical quantities). 
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Rules for calculating the mass operator and screened poten- 
tial can be found, for example, in [12] and [23] or in any other 
book on the theory of Green’s functions [24-27]. For our pur- 
poses it is sufficient to note that in the case of an interaction 
with impurities 


N 
U (r) = >; V (r—R4), (5.3) 


i=l 


where the index i enumerates the impurity atoms, R; is the ra- 
dius vector ofthe ith atom, andV (r — R.)is the screened potential 
of the interaction of an electron with the ith atom [for example, 
the Debye potential (4.2)]. 

The mass operator, generally speaking, is not Hermitian 
(which also indicates the presence of an attenuation of any ‘‘one- 
particle’? excitations). Therefore, together with Eq. (5.1) it 
would also be necessary to write a conjugate equation. For 
brevity we will simply bear this in mind. 

It is convenient to define the function JU, (r) as the periodic 
field of an ideal lattice with a vanishing small number of ex- 
citations (Such conditions are realized experimentally in an 
intrinsic or lightly doped semiconductor at sufficiently low 
temperatures, e.g., liquid helium temperature). 

Being interested in states with quite small crystal momenta 
(9><hb, where 6 is the reciprocal lattice constant) or quite 
shallow localized levels, we can generally eliminate U, from 
Eq. (5.1) with the aid of the well-known effective mass method 
[28-32]. In this situation the kinetic energy operatoris replaced 
by W,(— ih), where W,(p) is the dispersion relation in the band 
of interest to us, i.e., an eigenvalue of the auxiliary equation 


h2 


7 ?Xip + U Xp = Wi (P) Xo, (5.4) 


Mo 


pertaining to the Bloch function X,,.* Here mis the mass of the 
free electron. Thus, instead of Eq. (5.2) we have 


{E —Wi(— in) -U@}p—\dr' Mr, '; Eyer’) =0. 6.2’) 


From now on we will always utilize this equation (or its gen- 
eralization to the case of degenerate bands). 


*The extension to the case of degenerate bands is well known [28, 30-32] and is 
unimportant to us here. 
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If one approximates the function W,(p) by its expansion in a 
Taylor series around the band minimum and if it is assumed 
that the effective mass tensor is isotropic, then Eq. (5.2’) as- 
sumes the form 


(E+ 5 V?—U) ve) (dM (r, e's Eyp)=0, 6.2") 


where m is the effective mass. Replacing U(r) here by V (r)and 
neglecting the mass operator, we obtain Eq. (4.1). 

The first approximation—the replacement of the potential of 
all the impurity atoms by the potential of only one of them—has 
meaning if the screening radius is considerably smaller than 
the average distance between impurity atoms: 


rPyp<n Vs. 


This is nothing more than one of the conditions of light doping. 
The conditions for the applicability of the second approxi- 
mation—neglecting the mass operator—require a more detailed 
investigation. In the case of a degenerate electron gas the ap- 
propriate small parameter is [12] \)=e?/ehv,, where v, is the 
Fermi velocity. It can obviously be rewritten in the form 


A= (n2"?a9)}, (9.9) 


where ay = eh? /me? is the Bohr radius in the crystal. 

For a nondegenerate gas the parameter determining the rela- 
tive magnitude of the corrections associated with the mass 
operator is the usual Debye parameter, e?n}/*/eT (as usual, Tis 
the temperature in energy units) [33]. 

Let us note that neglecting the mass operator in the case 
given is not essential. Theoretically, taking it into account 
would change nothing in the considerations of Chap. 4. 

In the absence of impurity atoms and other structural defects 
(i.e., when U(r) =0) Eq. (5.4) allows one todetermine the energy 
spectrum of charge carriers in an ideal lattice while taking ac- 
count of the interaction between electrons. Aslongas the periodic 
field of the lattic is already eliminated thesystem being consid- 
ered becomes spatially homogeneous in the absence of the poten- 
tial U. The latter means in particular that M(r, r’; E) =M( r—r’, 
E). Let us represent the functions M and » in the form of 
Fourier expansions 
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M(r—r’,l4;E)= (dket —") M(k, l; E), 


f (r) = \ dke2™*" p(k), (5.6) 


where the wave vector k is related to the crystal momentum by 


(h is Planck’s constant, not Dirac’s). Then we easily find from 
Eq. (5.2) that 


E—W(p,)—M(k, f; E)=0. (5.7) 


This is the dispersion equation whose solutions determine the 
desired spectrum. Since the mass operator is not Hermitian, 
the function M (k, J, E) has both real and imaginary parts. Con- 
sequently, the solution of the dispersion equation has the form 


E=E, (k, |) — iE. (k, D, (5.7) 


where £, and -£o=y are real constants. 

Thus, states with given crystal momenta are, strictly 
speaking, not stationary. However, when 7% is small the imagi- 
nary part of M™ turns out to be small compared with the real 
part (this will be shown inChap. 7). This means that in a certain 
range of values of ne and T one can consider the states of the 
free charge carriers to be quasi-stationary and can speak of 
corrections to their energy arising from an interaction between 
electrons,* 

In the presence of a large amount of impurity, however, the 
Situation is changed. As already indicated in Chap, 3, the non- 
Stationarity of states with a given crystal momentum can prove 
to be important here, Consequently, the usual description of the 
system with the aid of a dispersion relation can become useless, 


*In addition to the rather trivial possibility considered here related to the small- 
ness of the interaction, there apparently exists still another one (34, 35]. That is, the 
imaginary part of the mass operator can be small, not “uniformly” in all momentum 
space but only close to the Fermi surface. Under the assumptionthat such is the case, 
this assertion can be proved exactly, independent of the interaction force between 
electrons [34]. On the other hand, the smallness of the attenuation justifies the as- 
sumption concerning the existence of the Fermi surface (35. 36], and in this sense the 
argument is self-consistent. However, it has apparently not been proved that states 
of the type considered represent a unique solution of the many-electron problem (even 
in the absense of the attractive force associated with phonon transfer). 
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A problem arises: what concepts and quantities should charac- 
terize the properties ofa system ofmany particles in the general 
case when there is no dispersion relation? 

An answer to this question can be obtained by formulating 
general equations for the macroscopic quantities of interest to 
us. Let us first consider the additive quantities whose average 
values differ from zero under equilibrium conditions, i.e., in 
the absence of external fields. Included, for example, are such 
quantities as the number of particles in the system, their kinetic 
energy, etc. 

According to Eq. (327°); in order to calculate such averages 
it is necessary to calculate the expression <d);a,,) . The latter 
can be considered as the limiting case of the time correlation 
function, <ap//(t’) Xdpr(t)), as t’->t. Here ap; (t’) and ap;(t) are op- 
erators in the Heisenberg representation; their equations of 
motion should be formulated with the complete Hamiltonian of 
the system. We will only be interested in equilibrium states. 
Then, as is well known, the density matrix has the form 


p = CexpB(uN — HH), (9.8) 


where 8=1/7, C is a normalization factor (determined from 
the condition Spur p= 1), Hand WN are the complete Hamiltonian 
and the operator of the total number of particles, respectively, 
and pv is the Fermi level. 

Under the stated conditions the correlation function depends 
only on the time difference, which is obviously physical and can 
be checked by a direct computation from Eq. (3.6) while taking 
Eq. (5.8) into account. Thus, one can write 


+00 
(pi (t’) Gpr (t)) = \ dye VENT (i, Is) (5.9) 


—oo 


(as usual, p = Ak). 

It can be proved that the spectral function /(k, /; v) defined 
by this relationship is real and nonnegative.* It is also seen 
that the amplitude A, Eq. (5.1), is 


+oc 
A= | dve tI (k, 1; vy eb (5.9") 


*We note that in [12] the spectral function I is introduced as the Fourier transform 
of the average value apt) apy (t)>. It differs from that used in Eq. (5.9) by the fac- 
tor exp B(E~—), where E = hv. 
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It follows from Eq. (5.9), in particular, thatthe energy spectrum 
of the charge carriers is determined by the poles of the spectral 
function. It is easy to express it in terms of some of the one- 
particle Green’s functions. Let us introduce, for example, 
retarded and advanced anticommutative Green’s functions: 


Ge(k, 1; £—t) = 0(E—F')< fa, apr tM» (5.10) 
Ga (k, ly f— t’) = — (0 (¢’ — t) ¢(api(2), Apl (0) Fox 


where 86 is the familiar step function 


Ly bot" 
t= |g op 

0, fr, 
and the averaging is carried out with the aid of the equilibrium 
density matrix, Eq. (5.8). Let us define the Fourier transforms 
of G,(k, /; v) and Ga(k, /; v) by equations of the type (5.9), 


-+co . 
Gralk,  t—t)= | dve@*">G, alk, bv), 


-—-00 


Then a simple calculation [taking Eq. (5.8) into account] shows 
that 


2/mG.(k, [; v) 2/mG, (k, & v) 
° —- ere ee FS _— L!:._._..\-_..._—— oe 
I(k, 5 v) FE 4 RED a (5.11) 


where E=Av, and the origin of the energy scale is chosen 
arbitrarily. 

In the absence of structural defects the Green’s function G,, 
is given by the expression 


hi 
Ad oe Ses 5,12 
in which the pole (if it lies on the real axis) should be bypassed, 


displacing E to the positive imaginary portion. Henceit follows, 
e.g., that for a real mass operator 


Im(G, (k, 1; v) = 5 o[hv — Wk, !)— M(k, ; E)) = 
A 6 fhv—Ex(, DI 


20, A OM (Kb ED’ (5.13) 
h dv 
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where, as before, the function £,(k, /) describes the dispersion 
relation calculated while taking account of the interaction be- 
tween electrons. 

According to Eq. (5.11), in the presence of the dispersion 
relation the spectral function also turns out to be a delta function. 

The spectral function is the first of the desired exact char- 
acteristics of the system. It is, however, somewhat awkward, 
because—even in the absence of the band index /, which can be 
dropped in the effective mass approximation and when the spin 
dependence of G,r,a is neglected—I depends on four variables, 
namely, the frequency and the three components of the wave 
vector. For many thermodynamic applications the density of 
states p(E) associated with this function proves to be more con- 
venient. In order to introduce it let us consider the expression 
for the total number of electrons N, inthe system. We make use 
of Eq. (3.7') by converting it from a summation over p to an 
integration and remembering that in the presentcase (p,/|A] p,)=I. 
We have 


Ne = 3) \ dp (Gap), 


rf 


where 2 is the volume of the system. Inserting here Eqs. (5.9) 
and (5.11), we obtain 


+co 
N 
ne =— = \ de—__°©_, (5.14) 
exp B(E — ph) +1 


where the ‘‘density of states,’’ p(E), is defined by 
9 (E) = = >) \ dkimG,(k, 1; v), E = hy. (5.15) 
l 


Equation (5.14) formally resembles a familiar equation from the 
statistical mechanics of gases with which there is also associ- 
ated a density of states function. For an ideal gas p(E)= po 
actually represents nothing more than the number of levels per 
unit volume and unit interval of energy. In this case we actually 
have M =0, the summation over / reduces to a summation over 
spin and gives a coefficient of 2, and, according to Eqs. (5.13) 
and (5.15), 


99 (E) = ~ \ dké (» = =). (5.16) 
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Introducing here the new independent variables W/h, a, and ), 
where a and b are some Suitable functions of k (for example, 
polar angles), we obtain 


po(E) = = \ dadbZ(a, 6, E/h), 5.16’) 


where Z (a, 6, W/h) is the Jacobian transform to the new varia- 
bles. This is the standard equation of the one-electron theory 
(see, for example, [3]). In particular, when W (k) = 4%2 /2m, the 
familiar equation 


a Gre 5.17 
Po (E) o— (2)2h3 0(E). ( ° ) 

follows from Eq. (5.16’). 
For a nonideal gas, but in the presence of a dispersion re- 


lation (ImM=0), the density of states, according to Eq. (5.13), 
has the form 


p(E) = 2\ dkd [E —W(k) — Mk, B)} (5.16) 


in analogy to Eq. (5.16). In the isotropic case (i.e., when W,and 
M depend only on |k|, we obtain from this 


OO 
a 
ro 


(5.17') 


where the function R(£) is the root of the equation EF — W (k)— 
M = 0. 

Finally, in the presence of a discrete level with energy E; 
the density of states is a delta function spike: in the region 
E~E;, 


0 (E) = N®(E — Ep: (5.17°) 


where JN; is the concentration of impurity atoms or other defects 
responsible for the formation of the given level. 

In the general case, when attenuation is present, it should 
be emphasized that the variable E in Eq. (5.15) is not the energy 
of an electron or any quasi-particle but rather just an integra- 
tion variable in Eq. (5.14). In this sense the term ‘‘density of 
states’? is not precise, since there can be no stationary states 
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with energy E. The function p(£), however, plays the role of the 
density of states in Eq. (5.14) andonly the following is necessary 
for obtaining the thermodynamic variables: after having calcu- 
lated p(£) from Eq. (5.15) wecan then relate the electron density 
with the Fermi level and the temperature by the aid of the 
standard equation (5.14). 

As seen from above, the density of states, as determined by 
Eq. (5.15), is an exact characteristic of the system. It can al- 
ways be introduced independently of whether a dispersion rela- 
tion exists and generally independently of what may have been 
assumed concerning the nature of the interaction between elec- 
trons. It is easy to extend Eq. (5.15) to the case of spatially 
homogeneous systems. In this case a Green’s function in the 
(p, /) representation is not diagonalized but, as seen from Eqs. 
(5.7) and (5.14), it should be in Eq. (5.15), its diagonal ele- 
ments being (1/2) G,(p, J; p, 4; v). In other words, Eq. (5.15) 
defines p(£)as the spur of the Green’s function (regarded as a 
matrix). This means that to compute p(E£), any suitable basis can 
be used. In particular, a coordinate representation can be 
utilized; then 


0(E, l) = — (dk Im G(k, t; ¥). (5.15') 


Sometimes it is convenient to introduce separate densities of 
states for the various bands. Obviously they are given simply 
by Eq. (5.15) in which the summation over / is omitted: 


2i ; W 
p(£) = + \ axl G(x, x; v). (5,15 ) 


A summation over spin is not performed in this equation. 

The concept of density of states is sometimes used for an 
exact many-electron definition of the concept of a ‘‘band.’’ That 
is, the /th band is designated as that region of energy in which 
the function (5.15”) differs from zero and is not singular. In the 
presence ofa dispersion relation this definition obviously coin- 
cides with the usual definition; however, in the presence of 
attenuation, the band, so defined, does not mean a region of al-~ 
lowed energies [the states(k, /) are not stationary]. 

The density of states, being an exact characteristic of the 
system, does not, however, yield all its properties. For ex- 
ample, as one may see from expression (3.7') for the average 
current density, in kinetic phenomena it is necessary to know 
the number of filled states n;(p, {) calculated under non- 
equilibrium conditions, i.e., in the presence of external fields; 
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on the other hand, the density of states is defined by Eqs. (5.10) 
and (5.15) as an ‘‘intrinsic’’ characteristic of the system in a 
state of thermodynamic equilibrium. 

The exact characteristics of the kinetic properties of the 
system can be determined by calculating the density matrix in 
the presence of any external perturbation and then formulating 
the average value as in Eq. (3.7). It is found possible to intro- 
duce a generalized electrical conductivity tensor that describes 
exactly the reaction of the system to a weak electric field [38- 
43]. That is, we let the field intensity have the form 


Esa iae (5.18) 


where Ema are constants, anda, B, are vectorindices. Then, to 
an approximation linear in £., the expression 


Gane (5.19) 
is valid for the average current density, where 
Visa = Oap (o) Em, Bs (5.20) 


and the electrical conductivity tensor ogg (w) is given by 


+00 
cap (o)= + \ 6O<Tie(®, dp OL > ede, (5.21) 


Here the averaging is accomplished with the aid of the equi- 
librium density matrix (5.8); d,(0) is the Heisenberg operator 
of the dipole moment of the system at time ‘=0. Integrating by 
parts, we may put Eq. (5.21) easily into the form 


-++00 
baat \ dt) <IjnO, mL De. (5.22) 


—0O 


° 
a 


Oap (w) = 


may 


Here m is a quantity with dimensions of mass (sometimes called 
the inertial effective mass) for which an explicit expression can 
be found (e.g., in[12]) (See also Chap. 12); ug (0) is the Heisenberg 
velocity operator for t=0. 

Sometimes a somewhat different form of Eq. (5.21) is con- 
venient, That is, by noting that 


Ad, (t) 
at 


= \ dx’ ig (x’, 2), 
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we obtain from Eq. (5.21) 
+00 i 
< 0 het , ° e , 
Gap (0) = lim \ dtO(t)e \ dx Cia (x, 4), ja(x’, OLD. (5.23) 


—0O 


In a spatially homogeneous system the correlation function of 
the currents, <[ja(X) 4), jp(x’, O)I_>, obviously depends only on 
x—x’ and ogg and is independent of the coordinates. 

According to Eqs. (5.21)-(5.23) and Eq. (3.7), the complex 
electrical conductivity is expressed in terms of a two-particle 
Green function of the form 


H (1) = 0) < laa, a(0)a(0)]_> (5.24) 


(for the sake of brevity, we have not written out arguments such 
as p, /). Generally speaking, the Fourier transform of this 
function may not be expressed only in terms of the density of 
states and is an independent exact characteristic of the system. 
In some cases, however, it is possible to establish an approxi- 
mate direct relationship between # and p (see Chap. 13). 


CHAPTER 6 


Equations for Green Functions 


As we saw in the previous section, the calculation of the 
density of states (and with it a number of thermodynamic prop- 
erties), as well as the calculation of the complex electrical 
conductivity, is reduced to determining Green’s functions in one 
form or another. Let us now consider the formulation of equa- 
tions for the Green functions. To this end, according to [12], 
one should differentiate the desired function with respect to the 
time ¢t, keeping in mind the familiar rule for differentiating 
operators: 


nee as © =[L(t), #].. (6.1) 


Here L is any operator, not explicity a function of time; # is 
the total Hamiltonian of the system. 

We will assume that the firststepin the solution of the many- 
electron problem—segregation of the conduction electrons and 
holes—has already been accomplished. Then the Hamiltonian of 
the system comprises the kinetic energy of the electrons (or 
holes), their potential energy in the periodic field, Uz, and in 
the field of the lattice defects, and also their Coulomb inter- 
action with one another: 


H= \ dx dx'a (x) or oa Aye +6 (x — x’) [Uc (x) + 
AU), (x) a(x’) + = \ dx dx'a(x)a(x')Ve(x — x’)a/x’)a(x). (6.2) 
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Here a(x) and a(x) are the Heisenberg operators in the configura- 
tion representation (until the averaging over the impurity co- 
ordinates is done this representation is more convenient than 
the momentum representation); Uo(x) is the potential energy in 
the field of the lattice defects unscreened by the free charges, 
V- is the Coulomb energy of the interaction of particles with 
each other (calculated, just as in the case of Uo, taking the 
static dielectric constant « of the lattice into account). The 
Hamiltonian (6.2) is written for a homopolar semiconductor. It 
is a trivial matter to extend it to the case in which free charges 
of both signs are present. 

We will consider, as an example, the retarded Green’s 
function 


G, (x, x’, (—t’)=i0 ((—t’) < [a(x, f), a(x’, £9), (6.3) 


(the function (5.10) is its own Fourier transform with respect to 
coordinates in a spatially homogeneous system). Employing 
Eqs. (6.1) and (6.2) we find 


{in i ve SUL) U4 (9) G,(x,x,f—)— 


= \dx"Ve (x — x”) ‘Ca (x”, x”, x, t,x’, (1) = — hb (x — x’) b(t’). 
(6.4) 


Here Go, is the second Green’s function: 


GoAX a XG GX, = 16 (t—t’) < [a(x”, Da(x”, Na(x, 0), 
a(x’, t’)],) (6.5) 


Differentiating G(x, x’; f—t’) with respect to ?¢’, we obtain in 
an analogous manner an equation similar to Eq. (6.4). As the 
boundary conditions for this equation one may take the usual 
periodic boundary conditions within a cube of volume 2. These 
conditions will generally not play a role when the specimen is 
sufficiently large and we are not interested in surface effects. 
The initial condition (i.e., the boundary condition with respect 
to the variable ¢f— ¢t’), according to Eq. (6.3), has the form 


G, (x, x’; {—t')=0 ({—t’ <Q). (6.6) 
The effect of the screening of the external field by the free 


charges, as well as all other effects of the interaction between 
electrons, iS contained in the last term on the left side of 
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Eq. (6.4). It is obviously convenient toisolate the terms respon- 
Sible for this effect by changing over from the unscreened 
potential U,)(x) to the screened potential U(x). To this end we 
note that in a system of particles with Coulombic interaction 
the potential energy of an electron in the screened field of the 
defects, U(x), satisfies the Poisson equation (in Gaussian units): 


VU (x) = — = (04x) Ape (a). (6.7) 


Here po, is the space-charge density of the structural defects. If 
the motion of the impurity atoms is neglected, po; is independent 
of U. 

The functional Ap;{U] is the variation of the charge density 
of all electrons except the given one in the field of the impurity 
atoms, An explicit expression for Ap, accurate to terms of any 
order of U can be derived by using (e.g.,) the method of [44]. 
We restrict ourselves to a linear approximation (which, in the 
classical region, corresponds to the Debye approximation). 
Then, according to [12], 


,, +90 = 
Ao, (x) = — r \ a0 (E 1) < [a(x Nalx, t), AHO) >. (6,8) 


The operator A # describes the interaction of electrons with the 
screened field of the impurity and has the form 


Ad = \ dx'U (x')a(x')a(x’). (6.9) 


Let us introduce the Green’s function 


H (x, x’, t—¢) = 0 (t —¢’)< a(x, Da(x, ), a(x’, ax’, Uy. 
(6.10) 


Substituting Eq. (6.9) into Eq. (6.8) and recalling Eq. (6.10), we 
obtain 


+00 
Ap, (x) = —\ dx’ \ aH (x, x’; )UX)= 
=< \ dx’o¢ (x, x't v) |ycol/ (x’). (6.11) 


Here X#(x, x’; v) is the Fourier transform of #(x, x’; ¢) formulated 
according to the method exemplified in Eq. (5.9). Thus Eq. (6.7) 
assumes the form 
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4ne? 


vu += ( dx’ # (x, x"; NUK’) = — on (x). (6.12) 


On the other hand, the unscreened potential Uo, clearly satisfies 
the equation 


Pi (x). (6.13) 


Let us introduce the Green’s functions for Eqs. (6.12) and (6.13), 
assuming, as usual ,* 


VxD (x, x’) + ae \ dx" 20, x”; 0) D(x”, x’) = — 6(x—x’), (6.14) 


eh 
VD q(x, x’) = — 8(x’ — x). (614°) 
Then, obviously, 
U (x) = — dx'D(x, x’) p¢(x’), (6.15) 
U, (x) = ae ( dx'D,(x, x’) 94 (x’), (6,15’) 


and the last term on the left sideof Eq. (6.4) can be represented 
in the form 


\ Go, (x”, x”, x, £5 X', LY Ve (x— x") dx" = 


4ne 


- St | dx {D (x, x”) — Dy (x, X”)} 09 (X") G(x, x3 E+ 


ne \ dx" dt"M, (x, x”; t—t")G, (x",.x’; t’—t’). (6.16) 


The mass operator M, (x, x”; ¢— t”), defined by this expression ,** 
describes the change in the dispersion relation and the attenua- 
tion of one-particle excitations because of the interaction between 
electrons. From Eq. (6.16) we obtain 


*The function D(x,x’) is nothing more than the nonrelativistic one-photon Green’s 
function, D(x,x,v) for v= 0{12]. The function K (x,x’, v) is clearly related to the 
polarization operator [23]. 

** This definition differs from that given in [23] in that here the terms responsible 
for the screening are explicitly separated out. 


EQUATIONS FOR GREEN FUNCTIONS 37 


M, (x, X’5 Vv) = \ dx" AX!’ Go, (x", x”, X5 X75 v) Ve (x — x”) X 
x GP (x, x’, vy) — 228K — x’) ( dx” {D(x, x”)— 
& 
— D, (x, x”)} pe (x). (6.17) 


for the Fourier transform of M, withrespect to time. The equa- 
tion for the Fourier transform of the Green’s function, accord- 
ing to Eqs. (6.4) and (6.16), assumes the form 


vx — U(x) — U(x)} G(x, x’; v) — 


h2 
{E+ 2M 
- ( dx" (x, x”; v) G(x", x'3 v) = — hid (x—x’); E=2nfiv. (6.18) 


(Here and subsequently we omit the index r since this cannot 
lead to ambiguity.) The boundary condition (6.6) is transformed 
into an analyticity condition for G,(x, x’; v)intheupper half-plane 
of the complex variable v. Analogously, the advanced function 
Ga (x, X*; v), Should be analytic in the lower half-plane. 

Obviously Eq. (5.2), which defines the energy spectrum of 
the system, is nothing more than a homogeneous equation as- 
sociated with an inhomogeneous problem (6.18). Therefore it 
could be said that the energy spectrum of the charge carriers 
is determined by the poles of the one-particle Green function. 
When attenuation is present, however, this statement requires 
an explanation. In the complex plane the Green’s function is 
usually defined by the equation 


420 | | 
Go) =s- J dv eh (6.19) 


=O 


(following from, for example, Eq. (5.11) and the analyticity 
properties mentioned above). 

It is seen that if-the spectral function is not singular (as 
happens when attenuation is present), the sole discontinuity of 
the right side of Eq. (6.19) is the branch line on the real axis, 
However, the delta function spectral functions are caused by 
poles on the real axis (a system without attenuation). The right 
side of Eq. (6.19) has no poles in the complex plane. 

The solution of this apparent paradox lies in the fact that 
the function G is multisheeted and the representation of Eq. 
(6.19) is valid only on one (the so-called ‘‘physical’’). Poles in 
the complex plane of interest to us lie ona ‘‘nonphysical”’’ 
sheet of the Green function. They can be formally defined, for 
example, as poles of the retarded Green function in the lower 
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half-plane. This means that the mass operator Mx, x’; v), ap- 
pearing in the dispersion equation (5.7) should be calculated by 
making the substitution 


Vv—> V—iE, (6.20) 


where e>0, «+0. 

The relationships (6.12)-(6.15) assume Somewhat simpler 
forms in the case when the system inthe absence of a screening 
charge is spatially homogeneous. Then, obviously, the function 
H (x,x; v) depends on only the difference x—x’, and can be 
written 


i (x, x's v= (dke*™ (Ke X=") ap (k, vy). (6,21) 


Introducing in a similar fashion the Fourier transforms of the 
potentials and the charge density, we obtain 


1 p, (k) _ 1 tk) 6.15. 
Oa ee U,(k) ee ( ) 
neh ’ ) 


Whence it is seen, in particular, that the magnitude of % (k, 0) 
when k=O determines the asymptotic behavior of the potential 
far away from the charges producing the field. That is, the 
quantity 


a = lim 2% (k, 0) (6.22) 


k-—+0 


is nothing more than the screening radius* [12, 23] (for example, 
the Debye radius, if the Debye approximation is valid). 

Equations (6.17) and (6,18) should be supplemented by the 
condition of system neutrality. It is clear that the latter has the 
form 


\ (x)dx =0, (6.23) 


where 


p(X) = Pr (X)-1 pe (X) 


*Later we will see that lim W (k, 0) <0. 
k —0 
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is the total charge density in the system (po, takes all electrons 
into account without exception). In a spatially homogeneous 
system, at which we arrive after exclusion of the periodic field 
and averaging over all impurity configurations, condition (6.23) 
can be transformed to a somewhat more convenient form, i.é., 


o (x)= — e(n—n,)--Ap (x), (6.24) 


where n and n- are the average concentrations of impurity atoms 
and electrons, respectively, and Ap(x) is the redistribution of 
the electron density caused by the total field of all the charges 
of the system.* Of course, in actuality Ap=0; it is sufficient, 
however, to impose the condition 


\ d xAp (x)=0. (6.23°) 


For Ap one can formally write down Eq. (6.11), where, how- 
ever, U(x’) now should represent the total potential Utotal pro- 
duced by all the charges of the system. 

Keeping in mind Eqs. (6.21) and (6.22) and noting that n<<, 
we obtain from Eqs. (6.23), (6.23') and (6.24) 


jZ=n (6.25) 
\d x’ Usotal (x’) = 9. (6.26) 


The latter relationship, taking Eq. (6.25) into account, can be 
rewritten in the form 


Ve (k) Ug (k) kao = | dx’ {Ve (x) +U, (x')} 
er Vege: (k)-+U (k) k=0 = \ dx’ {Ve, scr (x')+U (x' p20, (6.27) 
where V,, scr is the screened potential, whose Fourier transform 
is given by 


R2V _(k 
Ve, scr (k) as c| 


ae (6.28) 
neh Ut 0) 


Relationships (6.25) and (6.27) are the most convenient for the 
later formulation of the neutrality condition. 


*The difficulty with the self-energy in the nonrelativistic problem is easily elim- 
inated (using, e.g., the limiting A-process f109]). 
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Eliminating the periodic potential of the lattice with the aid 
of the effective mass method and assuming the band of interest 
to us (for example, the conduction band in an n-type material) 
is nondegenerate, we obtain, instead of Eq. (6.18), 


{E—W (— ihyx)-+-C—U (x)} G(x, x’, v) — 
— \ dx"M (x, x”; v) G(x", x’; v) =— hb (x — x’). (6.29) 


Here C is a constant replacing the periodic potential U, (x). Its 
choice is, of course, arbitrary; in particular, in Eq. (5.4) itis 
assumed that C=0. For later work it is convenient to retain C 
in an undefined form: we choose this constant such that the 
neutrality condition (6.27) is automatically satisfied.* 

Equation (6.29) is written in the standard ‘‘one-band’’ approxi- 
mation. It should be remembered that in the theory of heavily 
doped semiconductors this approximation can prove to be in- 
adequate. Actually, as we will see later, a calculation of the 
density of states near the middle of the forbidden band can be 
of great interest. In this region, obviously, it is necessary to 
consider the conduction and valence bands (we will call them 
first and second) on the same footing.** In other words, we have 
here ‘‘nearly degenerate’? bands which must be considered by 
the methods developed in [30-32]. We will assume that 
each of these bands, individually, is not degenerate and, more- 
over, that the width A of the forbidden band (in the absence of 
impurities) is considerably less than the distance from the in- 
trinsic extrema of the given bands to all other bands. Then we 
are dealing with a ‘‘two-band’’? model and the *‘smoothed’’ wave 
function of the effective mass method is atwo-row matrix: 


et 7 
LS 
Let us designate by W,(&) and W2(&) the functions which de- 


Scribe the dispersion relations in bands 1 and 2, considered to 
be nondegenerate, and let us designate the bottom of the first 


*What has been said, of course, does not mean that the periodic field of the lattice 
has any relationship to the neutrality condition. The use of the constant C for this 
purpose is nothing more than a convenient technical device. One would be just as 
successful by simply shifting the reference origin of the energy and introducing C as 
a “renormalization” constant. 

**The use of the term “valence band” does not mean, of course, that the standard 
“one-electron” meaning is attributed to it. Here it 1s simply the system of Bloch 
functions [solutions of Eq. (5.3)] used for the second quantization. 


EQUATIONS FOR GREEN FUNCTIONS 41 


band as the origin of the energy scale. Then by a standard 
method [32] we easily determine a set of equations for the 
auxiliary one-electron problem: 


(W, (kt) —E+U) pi +daphakptp,=0 (6.30a) 
(W2 (k)-E+U) Wo+dapkakatp =0. (6.30b) 


Here E are the energy eigenvalues, V is an operator describing 
the effect of external fields, and the term dag describes the 
‘‘interaction”’ between bands, being given by the equation 


ae he a 8B x. 1 : 
dap = — oe Di Pi hi ee) (6.31) 


The index / enumerates all the bands other than the first and 
Second, &: - the energy value at the bottom of the /th band, 
and p3, and p? are matrix elements of the operators of the a and 
B components of the momentum for transitions between the 
corresponding bands. According to what has been said above, 
lex] =>A. Therefore, when the spacing between the first and sec- 
ond bands is infinite (A — ©), the tensor dag becomes zero and 
the bands become independent of each other, as they should be. 

It is convenient to rewrite the set of equations (6.30) in 
matrix form, assuming 


WitW.=2W, Wi—W.=273, 
(dap + dag)Rakp = 211, (dap — dag) Rakp=2i72. (6.32) 


Forming half the difference and half the sum of Eqs. (6.30a) and 
(6.30b), we obtain 


(W —E+U) p+ yop=o. (6.33) 


Here, o is the Pauli matrix vector; y is a “‘vector’’ with com- 
ponents 71» Ye 13; and / is the unit matrix (later this will not be 
expressed explicitly). 

The Green’s function and mass operator in the present case 
are two-row matrices in terms of the ‘“‘band’’ indices 1, 2. 
Thus, instead of Eq. (6.3) we have 


Gy (x, x‘, ¢ — t’) = i (t—t’) < [az (x, B), a(x’, UY), i=1,2  (6.3') 


The equation for the Fourier transform, Gy;;(x, x’; v), has 
the form [see Eq, (6.29)]: 
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{E+-C—W (— ihy.) — U (x)} Gi; (x, x’, v) — (6.34) 
— 1 OjnGp; (x, X’5 Vv) — ( AX" Miz (x, X"; V) Ge; (x”, x’, VP = 
=— hd, ;6 (x = x’). 


When the ‘‘interaction’’ between bands is taken into account, 
the band number, understood literally, isnolongera good quan- 
tum number. Here, however, the secular equation (6.33) has a 
solution with two branches which are naturally called the first 
and second bands. When 730, and Y2>0, they transform adi- 
abatically into the ‘‘true’’ conduction and hole bands. 

The density of states in the first band is given, as formerly, 
by an equation of the type (5.15'): 


px (E) = = dk Im Gu (k, ¥) (6.35) 


In the second band, however, a ‘‘hole’? normalization of the 
energy proves to be more convenient. In order to make this 
changeover, let us note that the electron density in the second 
band is 


4 p2(E)[e? & —Y 4-1) dE = 
+00 
= { dEps (E) (1 — fe OP +1}. (6.36) 


The first term in the second line simply gives the total electron 
density in the second band and isofno interest to us; the second 
term, however, gives the holedensity, n,, the term that actually 
enters into semiconductor effects. 


Thus 
+00 
nn = ( dEpa(E)fe? "9 417, (6.37) 
where jth = —yu, and the density of hole states, p,, is 
2 
p, (E) = ra Goo (k, —v). (6.38) 


According to Chap. 4 just this function is of interest to applica- 
tions: quasi-stationary states can describe the holes in the 
valence band. 
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Let us now return to the equation for the function #4 (x, x‘, -—?’), 
necessary for calculating the screened potential. It is clearly a 
special case of the expression 


Hx, yi x’, y'3 t—’) = itt’) <[a (x, tha(y, t), a(x’, t’)X 
xaly’, UL», (6.39) 


which, according to Eqs. (5.21)-(5.23), is also necessary for 
calculating the complex electrical conductivity tensor, 

For later work it is sufficient to calculate the function (6.39) 
in the ‘‘one-band’’ effective mass approximation. Let us intro- 
duce the notation 


<a (x) a(z)a(z)a(y)|a(x’) aly’) S> = 
= i6 (t—t’) < [a(x, t) a(z, tha(z, thay, t), a(x’, t'ya(y’, t’) >, (6.40) 


function. We obtain u, using Eq. (6.1) and the Hamiltonian (6.2), 
just as before, by separating out the terms responsible for the 
screening effect explicitly: 


{E + W (— ihy.) — W(— ihy,) + U(x) -—U(y)} #& ys x ys v) — 
— {dz {Ve@—y)—Velz — x)} (a(x) a@a@aly) [ax )aly’)>>— 


— SE (ae {D (x, z) — D, (x, z) — Dy, z) + 


+ Do (Ys Z)} Po (Z) 0 (% Ys XT V)= 
= h {6 (x—y’) <a(x’) a(y)> — 0 (x’ — y)< a(x)aly’)}- (6.41) 


Then, according to the definition of Eq. (6.39), the desired func- 
tion of # should be analytical in the upper half-plane of the 
complex variable v. 

The last two terms on the left side of Eq. (6.41) describe the 
effects of electron interaction, which do not amount simply toa 
screening of the external field. 

Just as in the case of Eq. (6.4) [or Eq. (6.29)], Eq. (6.41) 
contains a Green’s function of higher order than desired. Differ- 
entiating it and continuing this procedure further, one can obtain 
an infinite network of equations, equivalent, in fact, to the 
well-known equations in functional derivatives (Schwinger equa- 
tions). For an explicit calculation of the energy spectrum, 
density of states, etc., it is necessary to solve this network by 
some—by necessity—approximate method.* This problem is the 
Subject of the following chapters. 


*In particular, at a temperature T = 0 one cannot calculate the associated retarded 
and causal Green's functions for which the usual Feynman graphical technique is valid. 


CHAPTER 7 


Perturbation Theory: Summary 
of Results 


In proceeding toward an approximate calculation of the 
Green’s functions for the case of a doped semiconductor it is 
convenient to consider first a simplified statement of the prob- 
lem in which the interaction of electrons, both with each other 
and with the impurity, is treated as a small perturbation, 
Actually, such a statement of the problem does not satisfy the 
requirement of Chapt. 4 because it is known that it is impossible 
in any finite order of perturbation to obtain the bound states if 
the spectrum of the unperturbed problem is continuous. And, 
in fact, we will see later that certain characteristics of a heavily 
doped semiconductor are impossible to obtain accurately by any 
expansion in terms of the impurity concentration or in terms of 
the ‘“‘smallness of the interaction’? of the electrons with the 
impurity. For this reason a calculation by perturbation theory 
can be applied more readily to the case of moderate doping when 
we are not interested in the impurity band but in the effect of 
the impurity on the structure of the conduction band. Neverthe- 
less, this calculation also hasacertain interest for the theory of 
heavily doped semiconductors. First of all, the electron inter- 
action, as we will see, can actually be considered by such a 
method. Second, a calculation along the lines of perturbation 
theory permits, to a certain extent, asubstantiation of the state- 
ments of Chap. 3 concerning the possible lack of a dispersion 
relation in a doped semiconductor. Third, and finally, it ap- 
pears that the equations derived below can be used for a calcu- 
lation of the Fermi level in a heavily doped semiconductor 
(Chap. 10). 
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The problem in which we are interested has been considered 
in [45-51]*—by somewhat different but essentially equiva- 
lent methods and with similar results. Here we will follow 
primarily [48, 49, 51] which used practically identical 
methods. 

The result of the calculation depends on the degree of de- 
generacy of the electron gas. We will consider the case of total 
or nearly total degeneracy by assuming 


>>, (7.1) 


where the Fermi level, u, is measured from the bottom of the 
conduction band.** Actually, in the case of moderate doping the 
Situation is somewhat different; a calculation carried out for a 
nondegenerate gas [48] shows, however, that the principal con- 
clusions reached below also apply there. 

When T = 0, the parameter that determines the smallness of 
the interaction between electrons is given by Eq. (5.5). The 
smallness of the interaction with the impurity is determined 
first of all by the condition of applicability of the Born approxi- 
mation and, second, by the familiar condition 


<i, (7.2) 


where o is the characteristic scattering cross section at an 
isolated impurity atom and k is the characteristic wave number 
of the scattered charge carrier. From dimensional considera- 
tions it is clear that when scattering occurs at charged donors, 
in the case of total degeneracy both these conditions are ulti- 
mately reduced to the inequality (already known to us) 


A=(a n'y 1<l, (7.3) 


which also characterizes the interaction between electrons. 
Actually the inequality (7.3) is quite poorly satisfied for most 
known semiconductors so that a calculation based onit has value 
only as a model. Of still greater significance is another cir- 
cumstance: even if the inequality (7.3) is satisfied, the criterion 


*A very similar problem was also solved in (52]. A somewhat different form of the 
perturbation theory was used in {113}, where the problem of the density of states in a 
heavily doped semiconductor was solved by the method of moments. 

**Here and in analogous cases later on the terms “bottom of the conduction band,” 
“top of the valence band,” etc., refer to an intrinsic material. 
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(7.2), satisfied for most electrons, is not satisfied for the slow- 
est of them. For this reason it is impossible to uSe the results 
of a perturbation calculation for studying the density of states 
near the band edge. 

Under the conditions of incomplete degeneracy the criterion 
(7.3) retains a meaning if the temperature is sufficiently low 
[inequality (7.1)]. The additional small parameter T/\: which ap- 
pears in this inequality does not alter the structure of the per- 
turbation series if it remains smaller than (4 (which we will 
assume below). 

Within the framework of the problem stated we are still re- 
stricted by the ‘‘one-hand’’ approximation of the effective mass 
method [Eq. (6.29)]. Since the conditions for the smallness of 
both perturbation operators, U and M, are determined by the 
same parameter A, in the first nonzero-order approximation 
one can calculate the mass operator while neglecting the effect 
of the impurity on the behavior of the electrons. In other words, 
it is necessary to substitute into the right side of Eq. (6.17) the 
Green’s function for free electrons [with the dispersion relation 
W (k)]. Thus, we arrive at the well-studied problem of a high- 
density electron gas. 

Accurate here to the higher orders of A the decoupling 


Go (x”, x", x} x’; v) cen. G(x, x’, v) — n(x", x) G(x", x’; ¥), (7.4) 
is valid, where 


n(x”, x) = <<a(x"a(x) >= \ dkeW 27 (K+ X"—¥) n (Ky) (7.5) 


are the elements of the first (‘‘one-particle’’) density matrix. 
The second term on the right side of Eq. (6.17) can be neg- 
lected in the present approximation. Thus we obtain 


M(x, y; v) = { dke™™ *—Y) M (k, v) (7.6) 


M (k, v) = neV.(k) \e=o — \ dk'V,.(k’) n (k’— k), (7.7) 


in which the function n(k) must be calculated for an ideal gas: 
this is simply the Fermi function 


n(k) = [exp a iz (7.8) 


where to is the Fermi level calculated for an ideal electron gas. 
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Fig. 4. 


The first term in Eq. (7.7) should not be taken into account, 
It vanishes by virtue of the neutrality condition (6.27). The 
second term describes the usual exchange correction, The 
diagrams,* presented in Figs. 4 and 5 (solid lines—electron, 
dashed lines—‘‘interaction lines’’) correspond graphically to 
the first and second terms in Eq. (7.7). A direct calculation 
[see below, Eq. (7.14)] indicates that when 7=0, the exchange 
term in Eq. (7.7) has a logarithmic singularity at k=kr, where 
Rr is the Fermi wave number. The origin of this singularity is 
related to the long-range nature of the Coulomb forces (the 
potential V.(k) has a discontinuity atk = 0). Fromthis it is clear 
that the singularity is removed when account is taken of the 
screening [53]. The latter is achieved by replacing (6.28) with 


RV. (k) 
LOS 
ke ——"_ K (k,v) 

hi 


ITE 


(7.9) 


which corresponds to the insertion of all polarization loops into 
the dashed line of Fig. 5. Thus, the mass operator assumes the 
form 


Wig = ( dk’n(k’ — h) 


neh (7.10) 


k/? — K (k’, v) 


As shown in [53], when 7=0 a better approximation can be 
obtained by approximating the function K(k’, v) by the constant 
_ jim K(k’, 0). Then, according to Eq. (6.21), 


M\(k, v) = — = \ dk’0(u, — W (k’ — k)) (7.11) 


4 
Ri? + (Qn)? 


*The correspondence should be understood in the sense stated in the footnote on 
p. 43: when T = O one can calculate the source Green’s function with the aid of the 
usual Feynman technique; in Figs. 4 and 5 are presented graphs which depict the 
contribution of the corresponding mass operators. [et us note that condition (6.27) 
excludes not only the graphs of Fig. 4 but also other graphs for which Fig. 4 is the 
skeleton. 
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ae oO on 


= ~ 


a A ee 
Fig. 5. 


Let us note that this expression |like Eq. (7.7)] is independent 
of v and is real. The latter, according to Chap. 5, means that in 
the assumed approximation attenuation is absent and one can 
speak of the energy spectrum of the charge carriers (calculated 
while taking account of the interaction between electrons).* That 
is, instead of the ‘‘raw’”’ dispersion relation W(k) the ‘‘re- 
normalized’’ relation emerges: 


W, (k) = £1 (k) = W(k) + M(k) (7.12) 


(we now omit the argument v in the mass operator). It is not 
necessary to take into account explicitly the term with the mass 
operator in Eq. (6.29).** 

In the isotropic approximation, when 


Wk) = —, (7.13) 


the integral in (7.11) is easily calculated. We obtain 


2e7k ke — ke + 2 ke +k)t+ x? 
A eR SS Popa A 

4kk pp (kp — hk)? + 2? 

k+kR k—k 
eiee (are tg a ae tg “I, (7.14) 
2k ya x 
where 
4 h2ki 

a 215 and y= 2m 


In accordance with our earlier statement, when x = 0 [which 
satisfies the approximation of Eq. (7.7)], there is a logarithmic 
discontinuity for k = kr. 


*It can be proved [12] that such is always the case in the first nonzero order per- 
turbation approximation if the interaction force is not retarded. 

** 4 pparently the same situation also occurs in a metal close to the Fermi surface 
(see footnote on p. 22). In a metal, however, the electron density 1s practically con- 
stant only because of the renormalized dispersion relation, and has a physical meaning 
(in static effects). On the other hand, in a semiconductor the electron density can vary 
within broad limits and both quantities, W and W,, have meaning. 
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When 7 +0, it is convenient to leave the mass operator of 
Eq. (7.7) in the form of a single integral 
1k + Rk’ 


(7,14) 


Let us now return to the perturbation produced by the im- 
purity. Introducing the renormalized dispersion relation, we 
can rewrite Eq. (6.29) in symbolic form 


(Gg°+U)G=1: (7.15) 


where 
Got = —W, —C —E) (7.16) 
2 


In accordance with Chaps. 3, 5 we have to investigate the 
zeros of the expression G-', where the symbol «~» means, 
just as in Chap, 3, an average over all impurity configurations. 


i ad 


~~ 


Obviously, G-!+G-'; therefore it is necessary first of all to 
determine either the function G by averaging it in accordance 
with Eq. (3.1), or to average Eq. (7.15) itself directly. Keeping 
in mind the application of perturbation theory, it is convenient 
to follow the second method. 

Let us first isolate the average potentials by assuming 


U=TU0+u', U=0, (7.17) 
and by transforming Eq. (7.15) into the form 
Go'G--U'G=1, (7.18) 
we have 
Gy = Gy" -- U. (7.19) 
Integrating (7.18), we obtain 
26,606, -GU GUC (7,18') 
Averaging this equation while taking Eq. (7.17) into account and 


including only the first nonzero terms of the expansion in terms 
of U’, we obtain: 
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lg =", 
G = Gy + GU'GU'G) + ..., (7.20) 
where the dots represent terms of higher order of smallness. 


Let us introduce the ‘‘mass operator’? M;, describing the 
effect of the impurity, by assuming 


Gt — Got eM;. (7.21) 
According to Eqs, (7.21) and (7.20), 


M; = —U'GU' +... (7.22) 


The operator M:z should not be combined with the Coulomb mass 
operator M (6.17). We note, however, that both of these op- 
erators are introduced for one and the same reason: they de- 
scribe, most conveniently, the variation of the functionG~! (i.e., 
ultimately, the energy spectrum) in the light of any interaction 
effects. 

In order to perform the averaging over the impurity con- 
figurations explicitly, it is necessary to define the probability 
density of the distribution of impurity atoms, P(Ri, .., Ry). 
appearing in Eq. (3.1). The exact determination is in itself an 
extremely complicated problem requiring a specialization of 
the model under consideration. It is clear, however, that ina 
macroscopically uniform medium the function P should be ap- 
proximately constant, if all the separations between impurity 
atoms |R;—R; |, (f#/j, i, j= 1, .., N) are sufficiently great. At 
small separations a correlation appears in the distribution of 
the impurity atoms. In particular, the incidence of twoim- 
purity atoms at the same point must be excluded, i.e., the 
function P must become zero when any two of the arguments 
are identical. In connection with this we note that the density of 
donors (or acceptors) in presently known semiconductors even 
at high doping levels is always small compared with the density 
of the normal lattice atoms (n < 10” to10?' cm-%), Therefore the 
probability of the incidence of two impurity atoms at one point 
is small even when correlation is neglected and it can be neg- 
lected by inserting suitable regions in the integral of type (3.1) 
(only if the average function is not excessively large for identi- 
cal or nearly identical arguments). In view of what has been 
said one can assume that a reasonable approximation to real 
life is achieved if the function P(R), ...,Rw)is simply approxi- 
mated by a constant.* Then, in accordance with the normaliza- 
tion condition (3.2), 


*With the aid of the functional technique an aberaging method is developed in [68] 
for an degree of correlation between the impurity atoms. 
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Pes ot: (7.23) 
oe tom 


In this case Eq. (7.22) in the configurational representation as- 
sumes the form 


Mi (x, x5 €—f) = — =F aR 2. ARyU' (x) U’ (x’) x 


y G(x, x’, £—0). (7.24) 


According to Eqs. (5.3) and (7.17), 


U’ (x) = > V’ (x — R,), 


f=1 


in which 


V(x —R) =VxE—R) — ay VaR _., dRy V(x —R) = 


— V(x —R)— —V (k) 


k=0 


Here V(k)is the Fourier transform of the screened potential of 
an individual impurity atom. Thus 


N 
U’ (x) = SV — R) — nV (kK) kao © (7.25) 


f=} 


Substituting this expression into Eq. (7.24) and performing the 
integration, we obtain (asymptotic when Nooc, Q+ 00, N/Q= 
noo); 


Mi(x, xX; t-—t)=—an \4 ke2ri(k. xx 1 V (kK)? G(x, x's £—1'). 
(7.26) 


Here it is assumed that V(—k) = V*(k). 

The Green’s function, averaged with the probability density 
(7.13), clearly depends only on the difference x—x’, Therefore, 
introducing the four-dimensional Fourier transforms of M, and 


— 


G in the usual manner, we obtain 
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M; (k, v) = —n \ dk’ \V (k’) |? G(k —k’; v). (7.27) 


Here we have omitted the average symbol on G because the un- 
averaged Green’s functions do not appear again in this section. 

According to Eq. (7.25), 

U = nV(k)| kao. 

In view of the neutrality condition (6.27) this term must be 
omitted (we recall that we have already discarded the corre- 
sponding term in the mass operator M). It is convenient to set 
the constant C equal to zero in this case. 

Substituting Eq. (7.27) into Eq. (7.21) and taking Eq. (7.16) 
into consideration, we derive an equation for the mass operator: 


LV (Ry)? 
E—W,(k— k’) — M (k —k’, v) 
(E = hv). 


Mi (k, v) =n \ dk’ (7.28) 


In the derivation of Eq. (7.27) essentially two approximations 
of different nature were made. First, it was assumed that the 
interaction of an electron with each individual impurity atom is 
small; second, a relatively low [in the sense of condition (7.21)] 
impurity density n was assumed.* The first of these limitations 
is easy to remove, Actually, terms of higher order in M; [sym- 
bolized by the dots on the right side of Eq. (7.22)] are of two 
types: one represents the correction to the scattering at a single 
center and is proportional to the first power ofn; the other 
allows for the interference of waves scattered by different cen- 
ters and is proportional to a higher power of n. Terms of the 
first type can be summed formally by introducing instead of 
|V(k)|? the exact scattering amplitude (see [54]). This sub- 
stitution can be helpful in taking account of short-range forces 
when, besides (5.5), other parameters can also appear. If, 
however, (as we will assume later on) V(£) simply represents 
the screened Coulomb field of a charged donor or acceptor, 
then the introduction of the exact scattering amplitude instead 
of the Born approximation is more meaningful: the difference 
between them ‘appears only when the inequality (7.3) breaks 
down, when the criterion (7.2) is also not satisfied. 


*Here and later in this section the dependence on n for a given electron density 
ne (on which the cross section o and the Fermi wave number also depend) is borne in 
mind. It is convenient in carrying out the calculation to assume these quantities are 
formally different and to keep condition (6.25) in mind only in the final result. 
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Within the assumed approximation framework Eq, (7.28) 
should be solved by integrations after having restricted the 
first of them. Keeping in mind the rule of Eq. (6.20), we obtain: 


= Mi = eee eee 7.29 
Re Mr (k, v= My =n \ dk =e , (7.29) 

Im M¢ (k, v) = Mz = an( dk’ | V(b’) PO[E — W, (k — k’) — 
— My (k —k’, v)}. (7.30) 


Substituting My; (k, v) in the dispersion equation (5.7), we 
find, under the assumed approximations, 


, Wk) = 
=~ ~ + = Wp, (k) 4- AV, 7.01 
Ey (k) = W, (k) tn \ dk TenWwEney =r AW, (7.31) 
M k, ? t ft 
t= ney = (ak |V (k’) 2 8 [WW (k) — W(k —k’)}. 
OAV pAv=W, (7 32) 


For an estimate of the screened potential produced by an 
isolated donor the usual Debye equation can be used.* Then 
under the assumed normalization 


e2 
= 7.33 
a se (k? + 2) : eee) 


Thus, in the isotropic case of Eq. (7.13) the integrals (7.30) and 
(7.31) are easily calculated. [In [50] the integrals (7.31) 
and (7.32) are calculated fer certain directions of k in the case 
of anisotropy of the type of n-Ge.] It is convenient to write sep- 
arate results for v>0 and v<0, because in the first case one can 
neglect the term M, under the delta function symbol in Eq. (7.30) 
and in the second case one cannot.** Let us introduce the sym- 
bols 


a, a=Ra, b=x0, Wr=— (7.34) 


Jag 


r—(Vb + Vo \/ nae tx —( Wb VP 
c= ASE yf ee vee 
(21)47 4 


*This is not obvious for a degenerate system ({55]), but is nevertheless possible 
if the doped semiconductor is under consideration (see Sect. 11). 
**We recall that the origin of the energy scale in Eq. (7.13) coincides with the bot- 
tom of the conduction band. 


PERTURBATION THEORY: SUMMARY OF RESULTS 55 


where, as before, 


_>, 38h 
0 mee 
Then we obtain 
a) when v > 0 
M, = 2% = (7.35a) 
; nm (x -—a-: 6)? + 4ab ° 
ae | a (7.36a) 
" mn? (x — a + 6)? 4- 4a6 
b) when v< 0 
M,=— “wy ; 7.35b 
ee We Vo A 
2na, c 76 (c) | 4na, { | 
M aoa ee ee eee fo ee (SS ° 
° m2 Ws (a + ¢ + 6)? —4ac (2n)4x% (c+(Vb + VY —x)?? 
(7 .36b) 


where 0(c), as before, is a step function. 

According to Eqs. (7.31), (7.32) and (7.35a), (7.36a) one ob- 
tains for the correction to the energy and the attenuation con- 
stant in the region v>0 the following equations: 


Wij Se ee es 7.37 
(“) xez WY 4 AW (ky ’ ( 
2nnet 1 

ee Wa ost 


where Wx= h?x2/2m is the energy of an electron with wave vector 
iP 

Something can be said concerning the correction tothe spec- 
trum in the literal sense only wheny<|AW/(&)|. In the given case 
this reduces to the inequality 


an ¥ 1. (7.39) 


The condition (7.39) can be satisfied only for very small k 
(where, incidentally, all calculations based on perturbation 
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theory become questionable). When k~k ~n “, the inequality 
(7.39) is completely invalid. On the other hand, under the con- 
ditions of the applicability of Eq. (7.3) it is reversed. This 
also means (Chap. 3) that the concept of a dispersion relation 
can be inapplicable if the corrections due to the interaction of 
charge carriers with the impurity are taken into consideration. 


SUMMARY OF THE RESULTS 
OF CHAPTER 7 


In the first nonzero approximation in terms of the small 
parameter A=(aon"s)", i.e., for a sufficiently high doping level 
and sufficiently high electron density, the interaction between 
electrons does not lead to an attenuation of states with a specific 
momentum. The total effect of the interaction between electrons 
can be described with the aid of the ‘‘renormalized’’ dispersion 
relation (7.12). In this instance the characteristics of the energy 
spectrum (for example, the effective mass) start to depend on 
the density of charge carriers; the form of the dispersion rela- 
tion itself also does not agree with the ‘‘bare’’ relation, which 
is valid for a very small number of electrons. 

Under the same approximation the interaction of charge 
carriers with impurity atoms leads both toa shift of their energy 
and also to an attenuation of states with certain crystal mo- 
menta. Near the Fermi energy the attenuation constant exceeds 
the energy shift appreciably. This means that when the inter- 
action with the impurity is taken into account, the concept of the 
dispersion relation of the charge carriers loses its meaning. 
Under such conditions it is impossible to speak of an effective 
mass in the literal sense of the word. If desired, one can speak 
only of ‘‘effective masses,’’ which enter into various directly 
observed quantities—for example, an ‘‘optical’’ effective mass 
(see Chap. 5), the effective mass for the density of states, etc. 
Since they are determined from different experiments, these 
quantities, generally speaking, do not agree with each other— 
even if they result from the isotropic model (7.13). 


CHAPTER 8 


Semiclassieal Approximation 


As noted in the foregoing section, the inequality (7.2) cannot 
be satisfied for slow electrons. The meaning of this inequality 
consists of the possibility of neglecting the interaction of an 
electron with several donors at one time. Takingthis effect into 
account probably involves the greatest difficulty in attempts to 
solve the problem with the aid of some form of perturbation 
theory (obviously an expansion in terms of the parameter no/k, 
the type used in the theory of multiple scattering, loses meaning 
when ns /k $1). 

In [56, 57] an attempt was made to get rid of this re- 
striction* by means of a numerical solution of the homogeneous 
Schrodinger problem with the potential energy in the form of 
an ensemble of randomly distributed delta functions. Ina heavily 
doped semiconductor, however, the delta~-function approxima- 
tion of the screened potential is not accurate. Actually it has 
meaning if the screening radius 7,=1/2nx is small compared 
with the average distance between impurity atoms, n7/%, With 
an increase in the electron density the screening radius, of 
course, decreases; however, in the absence of injection n-¥° 
also decreases. Under conditions of degeneracy of the latter 
the dependence is stronger. For an estimate of r, the familiar 
Debye equation for a degenerate gas can be used [59]: 


ro=3(5) V aon,-¥/8. (8.1) 


*See also the qualitative discussion in [58]. 


57 


58 ELECTRON THEORY OF SEMICONDUCTORS 


Whence, when n.=n, 
rons we (nays = a (8.2) 
0 D) 0 ) ’ e 


where 4 is the parameter of Eq. (5.5). The right side of Eq. (8.2) 
does not decrease but increases with anincreasein the impurity 
concentration and, consequently, under conditions of degeneracy 
the delta-function approximation is justified only for a suffici- 
ently low level of doping.* 

Following the work of [60, 61] we are restricted to a 
consideration of the asymptotic case of extremely heavy doping 
under conditions of degeneracy when the inequality (7.3) and 
the estimate (8.2) are valid. ** [In Chap. 11 it will be shown that 
taking account of the corrections to r, does not alter the in- 
equality (8.2) to an appreciable extent,] The Coulomb interaction 
between electrons can be taken into account by the method of 
Chap. 7, i.e., simply by replacing the ‘‘bare’’ dispersion relation 
with the ‘‘renormalized’’ relation (7.12). 

For the sake of simplicity, let us consider only the isotropic 
case 


2 


[This expression differs from Eq. (7.13) in that now the zero 
energy coincides with the ‘‘renormalized’’ bottom of the con- 
duction band.] An extension to the case of n-Ge type of aniso- 
tropy is not difficult but it is not necessary in the problems 
under consideration. 

Let us first consider the ‘‘one-band’’? approximation. Under 
the assumed approximations Eq. (6.29) and its conjugate equa- 
tion take the form 


{E +-C -1. — Vi — U(x) \Gq, x’, vy) = —hd(x—x’), (8.3a) 


{E i CO +} = gee O (x) G(x, x’, v)= —hd(x — x’). (8.3b) 


rf 
20 . 


t 
*In a compensated (but degenerate) sample, when n, <n, the product ro ns is clearly 
still greater than the right side of Eq. (8.2). On the other hand, in a nondegenerate 


gas when n, = n we have ro ns = (eT/4re2n4 JA, which can be markedly greater than 
unity 

** A similar approximation was used in [62, 108, 112]. The results of (112) are 
similar to those discussed above. 
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According to Eq. (8.2) it is impossible, under the conditions 
being considered, to isolate the ‘‘sphere of action’? of an indi- 
vidual impurity atom. At any point in the crystal the potential 
energy of an electron is produced by itsinteraction with several 
atoms at one time. For this reason the function U(x) changes 
slowly at distances of the order of n—"4, and the solution of Eqs. 
(8.3a, b) can be obtained by semiclassical means, i.e., with 
the aid of an expansion in terms of derivatives of increasing 
order with respect to U.* 

It is convenient to introduce the new variables 


R=*=* rax—x’ (8.4) 


Taking half the sum and differences of Eqs. (8.3a) and (8.3b) and 
neglecting higher derivatives of U, we obtain (a, B are vector 
indices): 


52 
(E+ ie+C tae Vi bag VR-UR)— 
4 0°U (R) =) ee 
— BTA RAR | HH RE) = — RO nae 
h2 e: aU (R 
{gs (wr 7g) = Fa mR, "| Ot R; E) == 0. (8.6) 


Heree >0,e—0. If the term ‘ce is introduced in the left side of 
Eq. (8.5), the following boundary condition is clearly indicated: 
the function G(r,R; £)is analytic in the upper £ half-plane, Equa- 
tion (8.6) can be considered as a supplementary condition im- 
posed on the solution of the basic equation (8.5). 

Let us choose a system of units in which hk = 2m, = | (the 
energy is measured in units of Ws = m,e+/2e2n?), Then the symbolic 
solution of (8.5) can be written in the form 


G(r,R; E) = i ( dse—est+é (E+C)s 141.8 (Pr), (8.7) 


U 


*Strictly speaking, U(x) is a slowly varying random function; therefore, considera- 
tions of the smallness of the derivatives have a literal meaning only after averaging 
with respect to the impurity configurations which will also be done at the proper time. 
In other words, the quantities discarded in Eqs. (8.5)-(8.6), etc., in themselves may not 
be small; they only represent a small contribution in the final equations - (8.20) and 
others. 


60 ELECTRON THEORY OF SEMICONDUCTORS 
where 


li = exp LS & TR —U (R)), 

a aU (R) 
ns expis( TP — Fels aR aR, | (8 .8) 
(in the assumed approximation thatv,commutes with 0U/OR,AR,) » 
The operators L, and Ll, are easy to ‘‘untangle’’ by making 


use of the smoothness of the function U(R). To this end let us 
consider the operator in a somewhat more general form 


L =expo(U + 7), (8.9) 


where ois a continuous parameter, U and 7 are certain opera- 
tors, and y-is a number, Let us set 


L = Sexp(on7T). (8.10) 


Then for S we obtain 


— =US + HIT, SL (8.11) 


with the ‘‘initial condition’? S(0)=1. For the commutator in 
Eq. (8.11) we have 


l 
[T, S]_ = evU+nr) >; 2 K,e-ont , (8.12) 
{>1 i! 
where 
KK, =(Aiy, U + nT), Ki = ([T,UJ_. (8.13) 
In our case T= — Vr, U=U(R). Then, neglecting derivatives of 


higher order than second, as well as the corresponding powers 
VU and V?U, we obtain from Eqs. (8.10)-(8.13) 


S = ew. exp{ —n(Sa+a)I, (8.14) 
where 


m=VU+ 2(7U,7)S—-Ai, 


0’U (R) a 


do = (VU) — 2 ae oR, aR,OR, (8.15) 
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In applying Eqs. (8.9) and (8.10) to the operators (8.8) it is nec- 
essary to set \ = —is, y = 1/4 in the first case and n = | in the 
second case. Furthermore, obviously, 


exp (7 8Vi Ae isv*) 8(r) = ( dk exp {2nikr +is(2x)2k?}. (8.16) 


It is easy to verify that the solution determined by Eas. (8.8), 
(8.7)-(8.10) and (8.14)-(8.16) satisfies the supplementary condi- 
tion (8.6). 

We obviously need to take intoaccount the derivatives of U(R) 
in the kinetics. On the other hand, in equilibrium problems (in 
particular, for the calculation of the density of states) one can 
be limited to the zero-order approximation with respect to the 
commutator, neglecting all derivatives of U(R) in Eq. (8.14). 
The relative error in the Fermi level associated with this is of 
the order of 4 for the neglected terms. The approximation ob- 
tained in this manner is clearly equivalent to the approximation 
Standard bands. 

In the approximation described we thus obtain 

G(r, R; E)=i \ ds( dk exp {— es + is[E+ C —U(R) — (20)°k*] + 
0 


+. Onikr}. (8.17) 


For the averaging of the Green’s function (8.7) using Eqs. 
(3.1) and (7.23) it is sufficient to compute the average from the 
exponential form 


N 
exp yi (R — R,), 


f=1 


where ! is an arbitrary function which goes to zero at infinity.* 
According to Eqs. (3.1) and (7.23) 


—— ‘ 
exp jR—R)=—- | Var’ exp f(R — R’) | = 
i=1 22 


= {! fe a dR fehRRD _ uy". 


*The terms in the exponent nonlinear in U readily reduce to linear terms, leading 
to a subsidiary integration. Thus, for example, 
—iss\? 


SY ls 
exp| (9 ,U) Jam \ a4 exp {—a4+ (=") (an -U)}. 
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In the limit as N -0oo, Q.,00, n'/Q =n< © this gives 


SS 
N 

exp 3} f(R—R,) = exp{n (dr [ef(r) — 1}. (8.18) 
i=1 


Applying the result of Eq. (8.18) to the function (8.17), we find: 


G(r, R; E)= \ dkeesare G(k, E); (8.19) 
where 
G(k, E) =1 {ds exp {— es + is[E + C — (2n)jk?] -+- a(s)} (8.20) 
and 


a(s)=n ( dr {exp[— isV (r)] — 1}. (8.21) 


According to the neutrality condition (6.27), the term containing 
\ drV(r) should be deleted from Eq. (8.21) [as it has already been 
deleted from the mass operator (7.7)]. From this it is clear 
that it is convenient to set 


C = — a’ (0). (8.22) 


Then the sum isC + « =a, in the exponent in Eq. (8.20) assumes 
the form 


ay (s)=n \dr {exp[— isV (r)] — 1 + isV (r)}. (8.23) 
We note that 
Rea, (s)<0, a1 (0) = a, (0) = 0, 
a, (0) = — 2a? = — 11 \ drV?2 (r) <0. (8.24) 


In particular, for the Debye potential, V= (Ze/er)e-r/re(Z is the 
donor charge in units of e), we have 


= ( anZ7e4 ro) (8.25) 


e2 
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When one utilizes Eqs. (8.20)-(8.21) in conjunction with the 
Debye potential, it is necessary to note the familiar condition: 
the semiclassical approximation is completely invalid if small 
values of r play a significant role. However, such a complica- 
tion does not arise here. 

In the ‘‘two-band’’ approximation the set of Eqs. (6.34) should 
be used instead of (6.29). As usual, under the conditions of (7.3) 
the Coulomb interaction can be taken into account by the simple 
substitution W+W,, and all the semiclassical calculations pro- 
ceed exactly the same as before. 

Let both bands considered be spherical: 


V,= usa W.= —A— 


2m, 2m, 


It is convenient to select a system of units in which A=2Q1)=1. 
Then, according to Eq. (6.82), 


W = — — + (Bi — Ba) (2a)? #2, 


where B,,. = ™,/m,,.. For the unaveraged Green function we obtain 
(under the approximation of spherical bands) 


G(r, R; E)= { dke2*G (k, R; E), 
G (k; R; E) =i \ dsexp {— es+is[E+C—W— isU(R)}- 


{cos 1s — iT sins }, (8.26) 
T 


where [see Eq. (6.382) ] 


r=VinP+inPe+y . (8.27) 


In averaging over the impurity configurations it should be 
considered that the two-band model can be of particular interest 
in the presence of compensation. Correspondingly, the averag- 
ing equation (8.18) must be extended to the case when there are 
several types of impurity atoms, Assuming, asusual, Eq. (7.28) 
is valid, it is easy to show that 


exp > f(R — Ri) = exp >) 7; \ dr {exp[— isV;(N}—1},  (8.18") 
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where the index j{ enumerates the various types of impurity 
atoms and 7; is the concentration of type-j/ atoms. 

In view of the neutrality condition (6.27), generalized to the 
case under consideration, it is convenient to choose the con- 
stant C in the form 


C=s hn, \ dr V;(r), (8.28) 
j 
and we finally obtain 


G(x, x’; E)= \ dker"r G(k, E), (8.29) 


G(k, E)=i \ dsexp {— es + is[E —W(k)] + a2 (s)}- 


0 


and 
{cos oi -— Sin rs} (8.30) 
if 
where 
O5/(S) == aI dr {exp [— isV,(r)] — 1 + isV;(r)}. (8.31) 
F 


We note, in particular, the case when the index / assumes two 
values, a and d (acceptors and donors), for which V,=—V¢g= 
V>0. Then Eq. (8.31) has the form 


Oy (Ss) = — (Ng |- Ng) \ dr {1 — cossV (r)} — 
— { (Ng — Na) \ dr {sV (r) — sin sV (r)} = tg + dary. (8.32) 


CHAPTER 9 


Density of States 


For computing the density of states in the ‘‘one-band’”’ approx- 
imation one can use Eq. (5.15), where by / is understood only 
the spin index. Since the function (8.20) is independent of spin, 
the summation over / simply results inafactor 2, and we obtain 
in the assumed system of units 


p(E) = = {dk Imi \ ds exp {— es + isE — is (2n)Pk* + a4 (5)}. (9.1) 


It is convenient here to interchange the sequence of the integra- 
tions with respect to k and s; however, caution should be exer- 
cised because the integral with respect to k diverges ass-0. 
We note, in this connection, the identity 


Imi {dk \ ds exp {— es -|- is (2)? k?} = 0, 


0 


Subtracting this identity term by term from Eq. (9.1), we obtain 
an expression in which the order of integration can be changed.* 
Thus we find 


t Se ds 
= ——_— op | == |. iEs — 9.2 
p(E)= a lme * p.f J sy Oxp {au (s) + Es—es}, (9.2) 


*Instead of this procedure one could introduce in the integral expression in Eq. 
(9.1) a “regularizing factor” of the type (-7k?), where 7 is a small positive quantity. 
The final result, of course, would not be altered. 
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where the symbol p./. designates the finite (as s,—-0) portion of 
the integral. 

In most equations p(£) actually enters under the integral 
Sign, and in calculating many directly measured quantities it 
would be much more convenient to leave p(£) in the form of Eq. 
(9.2). It is interesting, however, to determine the explicit form 
of the density of states in three characteristic regions of 
energy—near the Fermi level, near the bottom of the conduction 
band and deep in the forbidden band. 

When £>0, the integral in Eq. (9.2) can be calculated by the 
method of steepest descents. According to Eq. (8.29), a, Wsg= 
0(A75/4); we will also assume that pa (this will be justified by a 
later calculation—see Chap. 10), It is easy to see that the most 
important crossing point is s=0. A simple calculation gives (the 
final equations are written in the usual units): 


160 /2ah3 oe ee ae cae /4 : 7 
where 
gees (9.4) 
8a? 


and the /, are Bessel functions with imaginary argument, 
In the energy region E~p >a), Eq. (9.3) assumes the form 


e (E)=po(| Fs) (9.5) 


where f is the density of states in an ideal gas (p»>~ VE). Thus, 
the density of states is only slightly different here from the un- 
perturbed case. The reason is clear: with an increase in im- 
purity concentration the electron density also increases and 
consequently the screening effect is enhanced; as a result the 
average potential energy of the electrons in the field of the im- 
purity increases more slowly than the Fermi level. 

Near the bottom of the conduction band (i.e., when E~0) the 
density of states changes radically. According to Eq. (9.3), we 
obtain here 


3/e ak T FE 
E) (2mr) ** a '*PC/4) (1/4) 


xn! 2A3 


where I is a gamma function. 
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Equation (9.6) is valid for E<0 only if |E|<CE,=(me*/2e%n2) (nas), 
On the other hand, for E<0 and|£|>E£, the integration in Eq. (9.2) 
should be done by a more exact method if attractive forces be- 
tween electrons and impurity atoms are present, It is conveni- 
ent, first, to consider the model example by approximating the 
function V(r) by a spherical cavity or ‘‘hump,”’ i.e., by assuming 


V(r) = +80 (re —1), (9.7) 


where g and r, are positive parameters. 

Strictly speaking, it is impossible to apply the semiclassical 
untangling of the exponents in Eq, (8.7) to the combination of 
potentials of the form (9.7). It is clear, however, that the in- 
tegral which defines the function a,changes only slightly if the 
step function in Eq. (9.7) is slightly ‘‘smoothed’’ so that the 
summed potential varies slowly nevertheless. The quality ofthe 
approximation of Eq. (9.7) can be calculated in the familiar 
manner by noting that 


ais(s) =nvg(e*## —1 igs), (9.8) 


where v,= (4n/3)r?; on the other hand, for the Debye form of the 
screened potential the function 2,(s), oscillating, increases in 
amplitude, as soo, as In’s. The increase is associated with 
the singularity of V(r) at small spacings, which in the semi- 
classical solution of the problem should generally be cut off, 
Correspondingly the actual function «a,(s) should not differ too 
greatly from (9.8). 

In the approximation of Eq. (9.7) attenuation is absent and 
the spectrum of the charge carriers is given by the expression 


W = (2nyR? = nu og + Ig. 


Serious importance should not be attributed to this result, itis 
specifically a result of the nature of the approximation assumed 
and vanishes with a more precise treatment. Another essential: 
by substituting Eq. (9.7) into Eq. (9.2) and integrating with re- 
spect to s, we obtain 


e 


2 7 of : —_ 1 a 
ar >; a (E+ nog + Igy O(E + nug +!lg). (9.9) 
[==-0 


p(E) = 


If the interaction of an electron with an impurity atom exhibits 
the nature of an attractive force [the minus sign in Eq. (9.7)], 
then, as seen from Eq. (9.9), the density of states differs from 
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zero for all values of FE, including even E<0, As a matter of 
fact, for any E within the summation some value of / exists such 
that E—nuog+lg>0. In other words, as deep a ‘‘tail’’ of the 
density of states as desired is formed in the forbidden band.* 

The series on the right side of Eq. (9.9) is nothing more than 
an expansion of the density of states in powers of the impurity 
concentration. It is seen that it is impossible to obtain a “‘tail’’ 
in any finite series of a perturbation theory of this type (in ac- 
cordance, for example, with the results of [47]). 

It is easy to find the asymptotic form of p(£) when E<0O and 
|E|+co, To this end it is only necessary to note that under the 
stated conditions in the series (9.9) only terms with sufficiently 
large / differ from zero and, consequently, the factorial can be 
approximated by the Stirling formula and the sum replaced by 
an integral. In this way we obtain 


ln t's 
p(E)~ Vy exp HEI + y—(lm + Ye) In ~~} (9.10) 


nu» 
where J,, is a whole number just greater than |E|-+ nv,g, 


Cxy=ly —|E|—nyg <1. 


In the transition from the model example to the case witha 
more realistic potential, Eq. (9.9), of course, is no longer 
true. However, as usual the density of states does not become 
zero for any value of E if attractive forces are present between 
the electrons and impurity atoms. In the ‘“‘tail,’’ i.e., limit as 
|E| +00, we obtain 


= _ LEI | TEL 
op (E) ~ exp E, In "Eyif’ (9.11) 
where 
m_.e* mes 
f —'/; —— f 3 | é 
Es = 2 aT A =e are (n1ay)'/ ; (9.12) 


Equation (9.11) is formally valid if|E[S> Eo. 

It is easy to understand the reason for the formation of the 
‘‘tail.”? It lies in the fact that for a homogeneous distribution of 
impurity throughout the sample, there are localized fluctuations 


*In the presently assumed “one-band” approximation we actually do not go below, 
for example, the middle of the forbidden band. 
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in the number of impurity atoms in some selected volume. This 
leads to a localized variation in the electron energy compared 
with its average value, In particular, thereareregions in which 
the energy is decreased, which—in averaging over all configura- 
tions—appears in the form ofa ‘‘tail,’’ 

The fluctuational model provides an extremely graphic, 
though not rigorous, method for estimating the density of states 
in the tail.* Let us select some volume v,and we will assume, 
for simplicity, that the addition of each atom beyond the average 
number nvo, alters the electron energy by the same amount Eerr. 
Then the problem reduces toa calculation of the fluctuation prob- 
ability at which a given number, /, of impurity atoms appear in 
the volume v,.** The probability, denoted by P, is given by the 
Poisson equation 


P(l) = ad emt 
The density of states will be, obviously, 


p(E)~ SPQ VE + lEerr. 
L 


In the tail, fluctuations with /S> | are of special interest. This 
permits writing P in the form 


4 


P=-—— 
V 2x1 


exp | [— nvy — IIn - (9,13) 


Ny 


(making use of Stirling’s formula). 

Let uv, be the volume whose linear dimensions are of the or- 
der of the Debye radius r,. Then, clearly, Eosp ~E, = e2/ero and, 
moreover, nv,s>1 [in the light of Eqs. (7.3) and (8.1)]. In the 
region of energy not too far from the band edge small fluctua- 
tions exist when |/ —nv,|<nv,. The exponent in Eq. (9.13) can 
then be expanded in a series in terms of / — nup and we obtain 


(1 — nu)? ) 
Pw eee ot. 
exp | — 2NVp ~ Onv, f 


*Here we introduce a simplified variant of the arguments contained in (112]. 
**Here the density of screening charge carriers is assumed constant, nvg. In prac- 
tice the fluctuations in the number of impurity atoms involve a localized change in 
electron density that should reduce the density of states somewhat in the tail. This 
effect clearly becomes especially important if in its absence bound states can be 
created in the hole produced as a result of the fluctuation. We account for this very 
crudely by assuming that such bound states never appear. 
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According to the assumption 


({ — nv, = fai ; 


and, consequently, the density of states will have basically a 
Gaussian shape 


Heer) (9.11') 
FE? 


Cc 


where 
E, = E, V 2nv, ~ Wer *"* = Wa (na). 


This equation is valid regardless of the sign of the interaction 
potential (because the latter enters as the Square). In other 
words, it defines a tail near both bands—electron and hole— 
regardless of the type of doping impurity.* Its region of appli- 
cability, however, is restricted by the smallness condition for 
the fluctuations. Deeper in the forbidden band small fluctuations 
no longer play a role when the inequality If — nty|<<nvg is not 
satisfied. Either the volume uv, is appreciably smaller than the 
‘Debye’? volume, or the number / appreciably exceeds nr} [the 
latter is possible because, according to Eq. (7.3), the average 
Spacing between impurity atoms is n-3<ry], 

For /S=nv, and attractive forces Eq. (9.13) leads at once to 
Eq. (9.11). On the other hand, if repulsive forces exist between 
electrons and impurity atoms, then Eq. (9.11) is not obtained 
and the tail does not appear deepinthe forbidden band. Actually, 
in this case the electron energy is reduced with a decrease in 
the number of atoms compared with the average number, and 
this decrease cannot be less than n,. 

In the ‘‘two-band’’ approximation the situation is altered 
somewhat. That is, the densities of electron and hole states, 
Oo, and ?,, are given by Eqs. (6.35) and (6.38). The matrix ele- 
ments G,,; and Go appearing there are easily calculated with the 
aid of Eq. (8.30) by noting that 


Gi = — Sp (G -+ 63G), 


*Equation (9.11°) can also be derived from the general expression (9.2) by calcu- 
lating | for f’ <Q and approximating a,(s) by the first term in its series expansion 
in terms of s. 
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and 
| 
Goo = “9. Sp (G oe 03C). (9,14) 


Of most interest here is a study of the density of states a) ina 
compensated specimen and hb) when the ‘‘interaction’’ between 
bands is taken into account. 

Let us consider the firstof these problems in which the bands 
will be assumed to be independent (y,,y2—0). Substituting Eq. 
(8.30) in Eqs. (6.35) and (6.38) and performing the integration 
with respect tok, we obtain [Cf. Eq. (9.2)] 


pg(E)= ae -f- \ eS T9205) sein (E, -1-x5)—cos (Ey +o) }, 
1)? 3%? gs /? 
(9,14a) 
and 
OnE) = ori | Se 4 fein (Ex — 39) — 005 (Es 29) 
251) /? 5° s/? 


(9.14b) 


where a, and «x, are the real and imaginary parts of o., defined 
by Eq. (8.31). The origins of the energy scales in Eqs. (9.14a) 
and (9.14b) coincide with the bottom of the conduction band and 
the bottom of the hole band, respectively. 

Apart from the substitution 8: 8,, the difference between 
Eqs. (9.14a) and (9.14b) is the sign in front of x,, The meaning 
of this is clear: the charges of an electron and a hole are op- 
posite in sign (a, changes sign together with the potential V). 
Otherwise, Eqs. (9.14a) and (9.14b) are the same as Eq. (9.2)— 
as are all the results therefrom relating to the qualitative be- 
havior of the densities of states. Only twocircumstances should 
be kept in mind. 

First, now the resultant attraction or repulsion is deter- 
mined by the sign of the expression — (na — n,) V/r)(if this quantity 
is negative, then we have a resultant attraction for electrons; if 
positive, then for holes). As usual, the ‘‘tail’? is formed only 
in the case of a resultant attraction. 

Second, the neutrality condition (6.25) is now replaced by the 
equation 


Ng + Nh = Ng + Ne. 
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For this reason the density of free charges can be appreciably 
less than the impurity concentration, with a corresponding less- 
ening of the screening effect. 

The first of the factors mentioned clearly leads to a decrease 
in the density of states in the ‘‘tail’’? and the second to an in- 
crease in it. Which effect occurs depends on the specific shape 
of the screened potential and this will not be considered here, 

In changing over to take account ofthe ‘‘interaction’’ between 
bands, we note that the ratio 


Int irl | Fake dap |? 
a Eaaere, al oe 
3 ( tt | 
is small compared to unity only if the forbidden band gap is not 
excessively small. Thus, excluding the case of p-type Ge, one 
can consider the ‘‘interaction’’ between bands as a perturbation. 
Obviously it is sufficient to calculate p,/(E); the density of 
hole states is obtained from this by reversing the roles of donors 
and acceptors and by the substitutions E+—E, B:- foe. 
Let us introduce the quantity 


1 Z a 
4m p dQ | dapkakp [? =k, (9.1 D) 


where the integration is done over the whole solid angle. The 
quantity C/A plays the role of a perturbation parameter. In the 
first approximation with respect to ¢ we obtain 

Gis (ky v) = if de HEE) fet 


0 


— (2x)! je | seo oe — Cea mes |}. (9.16) 


The second term within the braces gives the correction to the 
case of noninteracting bands. To estimate its role it is suffi- 
cient to calculate the corresponding correction to the density of 
states in the ‘‘tail,’’ Ap-, within the framework of a model of the 
type (9.7): 


Va = — Va = — gO(r, — 1). 
We obtain, for E<0, |E|S>g; |A+E\|g: 


Ape == Ao’ J. Ap”, (9.17) 
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) C (0) 
Ap’ = ———..- — -gp,", 9.18 
p er gp (9.18) 
t ee (igt+ng) 
Apt = sb — A — Et B05 (na — rra)f/- 
Qn? (2nvon,)/*A2B,/2 
4 lo +1 9.19 
“exp {lp — (fo-+—} In og © ( 


Here p{” is the density of states for £ = 0 [given by Eq. (9.14a)], 


1, is a whole number slightly greater than 
A+ E| 
g 


Equations (9.17)-(9.19) are valid if Ap, < p!. 

The correction Ap’ does not play a significant role—it only 
reduces the density of states to some extent in the ‘‘tail.’’ The 
correction Ap”, however, is fundamentally significant. It differs 
from zero only in the presence of acceptors and, as is easy to 
verify, can lead to a non-zero density of electron states in the 
forbidden band even if there are no donors at all. This result is 
not too astonishing because the ‘‘interaction’’ between bands 
mixes electron and hole states (see Chap. 6, p. 40). It is never- 
theless of interest because it means that strictly speaking the 
‘‘tail’? of electron and hole statescan arise through the influence 
of any impurity—either donor or acceptor. It should be noted, 
however, that in a p-type specimen, for example, the develop- 
ment of an electron ‘‘tail’’ will bedetected simply as a decrease 
in the intensity of the hole ‘‘tail.’’ 

Thus, in a heavily doped semiconductor the concept of edges 
of the conduction and hole bands have, strictly speaking, a some- 
what conditional meaning. They are not simply displaced (be- 
cause of Coulomb interaction, for example), but are ‘‘smeared 
out’? because of the development of ‘‘tails.’* Thisis in essence 
nothing more than the emergence of impurity bands or, more 
precisely, that which was left of them as a result of taking 
screening into account. It is not extraneous to note that the 
screening is quite effective for ‘‘annihilating’’ impurity states. 
The function (9.11) decreases very rapidly with depthin the for- 
bidden band so that the density ofstatesis appreciable only near 
the bottom of the conduction band. Thetotal number of electrons 
in the region E<0 (i.e., below the bottom of the conduction band), 
n_, can be determined from an equation of the type (5.14): 


+- (Ng — Ng) Vo 


*In p-type Ge, apparently, the gap between the various valence bands is also 
“smeared out.” 
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0 


nt Meaney. 
‘ exp B(E— py) + 4 


Substituting here the density of states in the form of Eq. (9.2), 
we can (at T=0) integrate exactly with respect to E. The subse- 
quent integral with respect to s can becalculated approximately, 
neglecting the factorially small number of electrons in the ‘‘tail,’’ 
i.e., in the asymptotic region |E|S}E£). To this end, just as in 
the case (9.3), itis only worth approximating the function «;(s) 
by the first term in its Taylor’s series expansion about the 
point s=0,* 


a1 (s) == — as?, (9.20) 
Then it is easy to determine that 
no~ nh? =: n(nay (9.21) 


We also note that at the point E=—W, (the point where dis- 
crete impurity levels are located for light doping) the density of 
States (9.2) exists for no characteristics. 

An example of the experimentally verified consequences of 
the statement concerning the existence of a ‘‘tail’’ consists of 
the necessary appearance of the excess current in a tunnel 
diode [63]. As is known, the excess current has a whole series 
of causes: in particular, deep traps play an extremely important 
role [64]. When the latter are excluded, however, the excess 
current does not disappear completely. The portion that re- 
mains can apparently be considered as the normal tunnel cur- 
rent flowing ‘‘from tail to tail.’? A quantitative theory of this 
phenomenon is still lacking. Qualitatively this point of view is 
not inconsistent with experiments [65], according to which the 
minimum tunnel current point is shifted (for constant tempera- 
ture) to the right along the voltage axis if the concentration of 
the doping impurity is increased. 

Two comments should be made in connection with what has 
been said. 

First, from the foregoing [see especially Eq. (9.9)] it is clear 
that the development of a ‘‘tail’? is causednot by a strong attrac- 
tion of an electron to one isolated donor but by the combined 
effect of many impurity atoms located quite close to one another 
[this is possible when rp S=1~-'»; see (8.2)]. Consequently, just 


*lt is easy to see that this is equivalent to the perturbation theory of Chap. 7. 
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at this point an average with an uncorrelated probability density 
(7.23) can be unsound; correlations between impurity atoms un- 
doubtedly exclude the possibility of the development of very large 
aggregates. Thus the series (9.9) effectively terminates, to- 
gether with the ‘‘tail.’? One can imagine, however, that under 
conditions of heavy doping this circumstance is not very im- 
portant. A correlation of the type being considered obviously 
appears at distances of the order of the lattice constant d. Con- 
niga ed the number of atoms in an aggregate can amount, 
say, to (ro/d)*, and the energy to which the ‘‘tail’? penetrates 
into he forbidden band will then be of the order of (r,/a) Eo. 
Equation (9.11) will still apply within its limits if ros>d (actually 
for d= 4x107cem, d=5x10 §cm and x» = 10?! cm-3, r,/d ~ 4). 
It is seen from this that the density of states in this energy 
region becomes extremely small—comparable, apparently, to 
the density of states produced by any type of random structural 
defects such as dislocations, small-angle grain boundaries, etc. 
In other words the termination of the ‘‘tail’’ because of correla- 
tion under the conditions being considered apparently comes at 
the point where the function (9.2) has practically become zero, 

Second, although the conclusion concerning the existence ofa 
‘‘tail’? and its rapid decrease as |E'— « apparently does not give 
rise to doubts, the question ofits exact asymptotic shape is nev- 
ertheless ultimately unclear. The point is that in calculating 
o(£) in the region of large negative energies the assumed semi- 
classical approximation becomes sufficiently convincing, For 
this reason it is desirable to perform an experimental study of 
the density of states in the forbidden band. One way to do this 
could be an investigation of the excess current in tunnel diodes; 
another possibility will be discussed in Chap. 13. 


SUMMARY OF THE RESULTS 
OF CHAPTER 9 


As a result of the interaction of charge carriers with the 
impurity atoms the density of states in a heavily doped semi- 
conductor does not become zero at the bottom of the conduction 
band (or at the bottom of the hole band)—the region of ‘‘allowed’’ 
energies is shifted downward. Sharp band edges disappear and 
a ‘‘tail’’? in the density of states, extendingdirectly into the for- 
bidden band, appears. The latter concept thus acquires a some- 
what conditional meaning butnevertheless can be retained because 
the density of states in the ‘“‘tail’? decreases rapidly (faster than 


76 ELECTRON THEORY OF SEMICONDUCTORS 


exponentially) with depth into the forbidden band. The charac- 
teristic energy Fy), which determines the rate of decrease, is 
given by Eq. (9.12). As a result of screening the number of 
states in the ‘‘forbidden’’ band turns out to be relatively small. 
According to Eq. (9.21) itis apparently less than the number of 
impurity atoms. Deep in the conduction band (in particular, 
near the Fermi level in an uncompensated specimen) the density 
of states differs only slightly from the unperturbed value. 

A ‘‘tail’?? is formed, strictly speaking, in the presence of 
any type of impurity—both attracting and repelling electrons. 
The former, however, produces a greater effect, 


CHAPTER 10 


The Fermi Level in a Heavily 


Doped Semiconductor 


In the general case the Fermi level is defined by Eq. (5.14) 
in which the electron density must be calculated with the aid of 
the neutrality condition, For an uncompensated specimen, neg- 
lecting the intrinsic conductivity, n, is equal to n and conse- 
quently Eq. (5.14) assumes the form 


+00 
n= (dE. ceereee 2 (2 eee (10.1) 
Bs exp 3(E —pw)+4 1 


On the basis of the results of the preceding section one can 
imagine that the Fermi level will differ slightly from the un- 
perturbed value jto calculated simply for an ideal Fermi gas of 
density n with the dispersion relation W (k). Under conditions of 
smallness of the parameter 2? the Coulomb interaction leads to 
only a small renormalization of the energy and the interaction 
of electrons with the impurity alters the density of states notice- 
ably only where it is small. Thus one can write 


= Uy + Apayat+ Apimp » (10.2) 


where Attcou. and Atimp are the corrections caused by the effect 
of the interaction between electrons and the interaction of elec- 
trons with the impurity, respectively. In the approximation as- 
sumed these can be calculated independently of each other. 

In order to determine Ap.ou for T=0O we write Eq. (10.1) 
twice—once for the system being considered and once for a 
hypothetical gas of noninteracting particles (of the same density). 
We obtain 


77 
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ji Oe) 

\ dEp(E) = \ dEp, (E). 

Keeping in mind Eqs, (7.12) and (5.16) and restricting ourselves 
to an isotropic system, it is easy to determine from this that 


Atleoul = M (A) |e=ep- (10.3) 


where ke = (3n/8n) is the Fermi wave number.* Substituting 
here Eq. (7.14) for the mass operator and noting that in view of 
Eqs. (7.8) and (8.1) x < kr, we find 


2 


2e7k 
Abeoul = — ——— = 0(A). (10.4) 


& 


In calculating the impurity correction it is convenientto sub- 
stitute into Eg. (10.1) the density of states in the form (9.2). 
Then, for 7T =O one can integrate exactly with respect to — and 
the relationship (10.1) assumes the form 


on'/sh3 ut 


*/2 4 ran? 1 
n=) im] — temp f Cr anelii 1 (40.5) 


Here p’=p,+Apinps the function «, (uw?) is obtained from a, (s) by 
replacing s by uw?. The integral in (10.5) can be calculated with 
the aid of the approximation (9.20). Keeping in mind also that 
according to Eq. (8.25) and (8.1) a/p =0(A), we find after a 
simple but somewhat lengthy calculation 


= (2mp')/* 3a? 
= Pay” (14. 2). (10.6) 


It follows from this that \u,,, < 0 and 


AMMMimp 0 m *) 


Ito 
i.€., aS A approaches zero asymptotically the Coulomb correc- 
tion to the Fermi level plays a major role. It is easy to under- 
stand the reason for this with the aid of the results of Chap. 9. 


That is, the interaction of the electrons with the impurity affects 
the location of the Fermi level for two reasons: first, it alters 


*iquation (10.3) is a special case of the general theory of [66]. 


THE FERMI LEVEL IN A HEAVILY DOPEO SEMICONDUCTOR 79 


the density of states within the conduction band (near the Fermi 
level itself); Second, it leads to the appearance of a density of 
states below the bottom of the band. The first effect, according 
to Eq. (9.5), gives a correction of the order of a?/y’? = 0(A3); the 
second, according to Eq. (9.21), gives just the correction in 
Eq. (10.6). 

The temperature corrections (for w,>T) can be calculated 
in an analogous manner. Taking Eq. (7.14) into consideration 
and restricting ourselves only to terms of lower order in 4, we 
finally obtain 


ae 12 na.) a 
=n, {1—0.2 nay — 2. E [14 BO | t, 
6 | 0.2 (na5) Oo ae l =? (an) - 
(10.7) 
where 
Z (3m2)/*h2n°/s 
i a 2m ° 


Thus, for T<jto the thermodynamic properties of a heavily 
doped semiconductor with homopolar conductivity are almost the 
same as for an ideal degenerate Fermi gas. Let us note, never- 
theless, one characteristic difference of the given system from 
degenerate semiconductors with moderate or light doping. As 
is known (see, for example, [67]), in the latter case two 
degeneracy temperatures occur, an ‘‘upper’’ and ‘‘lower.’’ 
The appearance of the latter is caused by the ‘‘sweeping out’’ 
of carriers from the conduction band with their transition into 
the impurity levels. In a heavily doped semiconductor the low 
temperature degeneracy is absent bacause, as a result of the 
screening, the density of states in the forbidden band is negligi- 
ble. The upper degeneracy temperature 7 *, can be estimated 
roughly from Eq. (10.7), it being defined as the temperature at 
which the Fermi level coincides with the bottom of the conduc- 
tion band in an intrinsic material. Under these conditions, of 
course, the inequalities used before are no longer valid and 
Eq. (10.7) can only give an order of magnitude. In this way we 
obtain 


T* ~T5(1 — 0.23). (10.8) 


Here T )= 12 n*/yo is the degeneracy temperature calculated neg- 
lecting the interaction of electrons with the impurity, 

In a compensated specimen, containing equal numbers of 
donors and acceptors, the Fermi level is of course situated 
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within the forbidden band. Its exactcalculation hardly makes any 
sense at the present time because it is highly dependent on the 
shape of the density of states in the ‘‘tail.’? One circumstance, 
however, deserves attention. Because of the presence of elec- 
tron and hole ‘‘tails’’ in a heavily doped semiconductor—even 
under conditions of compensation and at low temperatures— 
there can be electrons (or holes) located above (below) the Fermi 
level.* Their density is small (factorially small), but different 
from zero nevertheless. One might suspect that this would lead 
to a finite metallic-type conductivity; however, the question 
concerning the participation of the charge carriers in the ‘‘tail’’ 
in transfer processes requires a Special study. 


SUMMARY OF THE RESULTS 
OF CHAPTER 10 


In an uncompensated heavily doped semiconductor the Fermi 
level is Situated within the conduction band. This assertion is 
valid when intrinsic conductivity is neglected, i.e., under condi- 
tions of heavy doping, at practically any temperature. At suffi- 
ciently high impurity concentration [(naj)*S> 1] and sufficiently 
low temperature (J<wyp, where the Fermi level is measured 
from the bottom of the conduction band in an intrinsic material) 
the Fermi level differs slightly from |, the Fermi level in an 
ideal degenerate gas [see Eq. (10.7)]. Unlike lightly and mod- 
erately doped specimens, a heavily doped semiconductor has 
only one (the upper) degeneracy temperature; it differs slightly 
from the degeneracy temperature of an ideal gas. 

In a compensated specimen at any temperature down to7=0 
there are electrons or holes located, respectively, above or 
below the Fermi level, though in small numbers. 


ok ; é é 
The exclusion comprises only the case of exact compensation: ng = ng. 


CHAPTER 11 


The Momentum Distribution Function 


For 7T=0 and a Screened Potential 


The statistical characteristics of asystem of charge carriers 
are certainly not described merely by one Fermi level. The 
momentum distribution function of the charge carriers can also 
be of considerable interest. Here it is sufficient to limit our- 
selves to the ‘‘one-band’’? approximation of the effective mass 
method (if the case of the degenerate bands of p-type Ge is ex- 
cluded). Then, according to [12]* 


Im G,(k, E) 


_ i 
nk)=— | dE Eat 


—_ 


> p=hk (11.1) 


where G, is, as usual, the retarded anticommutative Green’s 
function. In the problem with an impurity theimportant thing is, 
of course, the Green’s function averaged over the impurity 
configurations. 

If there is no attenuation and the mass operator is independ- 
ent of E, then, according to Eqs. (5.13) and (11.1), a Fermi 
function is obtained for n (k) with renormalized energies. How- 
ever, in the presence of attenuation the location of the Fermi 
level is altered and the function n (k).no longer has a universal 
form. Even at T=0 the Fermi ‘‘step’’ can be smeared out. 

Substituting the Green’s function (8.20) into Eq. (11.1), as- 
suming T=0 and integrating with respect to E, we obtain 


*E quation (11.1) differs from that given in [12] by a constant factor. This is asso 
ciated with differing definitions of the Fourier transforms in (12] and in this paper. 


8 1 
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Bley sis Im ( ds .—estis(u—W, )-+ar(s) (11.2) 
2 n Ss 
0 

The integral appearing here has no Singularity at s=0, 
Ima,;=0 when s=0. 

Sometimes another representation for the distribution func- 
tion is convenient. That is, let us introduce the total mass op- 
erator M which describes the interaction ofcharge carriers with 
each other as well as with the impurity (and, possibly, with 
phonons, etc.). The general rule for formulating such a mass 
operator can be given by taking the average (as in Chap, 7) over 
all the impurity configurations just in the equations for the 
Green’s functions. It is convenient to rewrite the infinite set of 
correlated equations thus obtained in the form of an equation in 
functional derivatives [68]. Then a formally closed expression 
for M is obtained which then permits a step-by-step calculation 
with the aid of one method of approximation or another. In par- 
ticular, under the approximations of Chap. 7 there is obtained 


~~ 


M=™M:-+M. 
By definition, the inverse Green’s function is 
G4(k, E)= Gy*(k, E)+ Mk, E), (11.3) 
where | 
Go'=Wi(k)—E. 


Substituting Eq. (11.3) into Eq. (11.1), we obtain (for 7=0)* 
<a. 
x.  (E—W—Re M+ (mM) 


For sufficiently low attenuation, M(k, £) can bereplaced here 
by M(k, W(k). Then [48] 


Im M (k, E) (11.47) 


sl _¥-V®) 11.4 
n(k) eee arctan, reat (11.4) 


where 
W (k) = W(k) + ReM(k, W(k) 


*It can be proved in a general form [69] that the imaginary part of the mass operator, 
in contrast to the retarded Green's function, is not positive. 
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is (approximately) the dispersion relation calculated taking ac- 
count of all interactions. 

Equations (11.4) and (11.2) can have the usual range of appli- 
cability: when ’<1, this is the allowed band remote from its 
edges. In this case ImM = Im M; because the imaginary part of 
the Coulomb mass operator M is a quantity of higher order in A 
(Chap. 7). Comparing the right sides of Eqs. (11.4) and (11.2), 
we obtain* 


Him M,(k, W (@)| =——=—#= a.) 


tan Im \ GS exp {—es+t is(iu —W)+ a,} 
0 S 

Equation (11.5) permits a judgment to be made concerning 
the attenuation of states with the crystal momenta p=/k. In the 
approximation assumed the corresponding attenuation constant 7 
is simply 1/h |ImM,(k, W (k))|.** 

The explicit calculation of the integral with respect to s ap- 
pearing in Eqs, (11.2) and (11.5) is difficult. We note, however, 
that consideration of the corrections to the spherical band ap- 
proximation leads, according to Eq. (8.14), to an effective 
termination of the integral at large values of s. Therefore, we 
may rightly approximate the function a, (s)by Eq. (9.20). *** Then 
we easily obtain 


| uw — W (rR) 
=— J] t on 11.6 
n(k) =~ {1+ Ef Hh (11.6) 
and 
[ue —W, (&) | 
Im M; (k, W (k)) | = 7, 
AW f-Bry eS 
2 2a 
where 
Erfx = \ e—¥* dy. 
a 0 


*For an independent determination of the mass operator from the rule in [68] Eq. 
(10.5) is, of course, converted to an identity. 


**To avoid ambiguities we note that the time +! is not in general the same as the 


relaxation transport time which determines the charge-carrier mobility (Chap. 12). 
***This is nothing more than the evaluation of the integral by the method of steepest 
descents, which is valid in view of its convergence and the large value of a, Eq. 
(8.25). The principal crossing point is s = 0. 
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In particular, whenW, (2) =yu(in the absence of doping this would 
be the equation defining the Fermi surface), we have 


Za 


hZ 1 le 4 
7 : V+ cy nile. (11.8) 
AVn uma, 3 


Here Eq. (8.25), in which the unperturbed value of (8.1) is sub- 
stituted for rp, was used. The last expression in Eq. (11.8), un- 
like Eqs, (11.2) and (11.5), does not allowthe limit n=0 because 
it is derived under the assumptions of ~<1 and n=n,.. Sucha 
limiting process should be formally performed in Eqs. (11.6), 
(11.7) and (11.8) by setting nto zero for a given screening radius 
ro (i.e., for a given electron density). 

It is seen thatasaresultofthe interaction of charge carriers 
with the impurity the Fermi distribution is ‘smeared out,’’ even 
at absolute zero. Inother words, the representation of the Fermi 
surface in a heavily doped semiconductor has no unambiguous 
meaning.* This circumstance should not be surprising because 
(Chap. 7) in the presence of doping the assumed dispersion rela- 
tion becomes inexact. 

The ‘‘smearing’’ effect, described by Eq. (11.6), is small 
under the conditions being considered: according to Eqs. (8.25) 
and (8.1) the second term in Eq. (10.6) rapidly becomes 


1V-=p)= 


lp— W, (0) 


with the removal of Y, from the Fermi level. This effect, how- 
ever, is significant in calculating the screened potential, to 
which we now turn. 

According to Eq. (6.15), to calculate the screened potential 
in a Spatially homogeneous system it is necessary to find the 
Fourier transform of the two-particle Green’s function (6.10), 
which is a special case of the function # (x, y; x’, y’; v), Eq. 
(6.39), 

By calculating # when 4 is small, the Coulomb interaction 
between charge carriers can be taken into account with the aid 


*This statement apparently contradicts the familiar theorem of Migdal 36]. In the 
present case, however, the conditions under which it is developed are not satisfied. 
Thus, for example, the familiar statistical weight argument that demonstrates the dis- 
appearance of attenuation at the Fermi surface, fails: the impurity atoms are treated 
as infinitely heavy “field sources,” so that scattering at them is completely elastic. 
The accuracy of this latter approximation is determined by the parameter m/M, where 
M is the mass of the impurity atom; within the framework of the approximate theory 
developed here this approximation is necessary. 
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of an approximation method—the replacement of W(k) by W,(k)— 


already familiar to us. Then, the solution of Eq. (6.41) can be 
represented in the form 


K(x, yi x’, Yo v= {dz {g (x, y; y'3 Z; v)n(x’, z)— 
— g(x, y; Zz, x; v)a(z, y’)}, (11.9) 


where, as before, 
n(x, Z) = (a(x) a(z)), 
and the function g satisfies the equation 


h2 
2m, 


{ hv + ie — F-vEt ga-Vd TU ®) — UY) bece yix' 9) = 


= hd (x —x')6(y—y’) (11.10) 


and the equation conjugated with it. 

The representation (11.9) is convenient because Eq, (11.10) 
has the same structure as Eq. (8.3) and can be solved by a sim- 
ilar method. In particular, in the spherical band approximation 
the potential energy U generally falls out of the solution. The 


effect of the impurity on the function # appears in this case only 
in the form of the elements of the one-particle density matrix 
n(x, z)* These must be calculated independently, in terms of 
the one-particle Green’s function G (z, x; v), from formulas like 


Eqs. (5.9) and (5.11) [12]. We note that to calculate % it is Ob- 
viously necessary always to average the product gn over the 
impurity configurations and not each of the factors separately. 


Therefore, in general, the function # is not expressed directly 
in terms of n(k) The exception is only the case of spherical 
bands (satisfactory for the screening problem), Within the frame- 
work of this approximation we obtain, on the basis of Eqs. (11.9), 
(11.10) and (7.5), 


K(x, x's vy =lim# (x, y; x, ys v= \ desi Xx") (kK, Vv), 
x 
y’—>x’ 


(11.11) 


*From this it is especially easy to see that spherical band approximations are in- 
adequate in kinetics problems. 
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where 


n(k’)— n (k’ — k) 


#(k, Vv) = aa see cnr COM ae 


The factor 2 enters here as a result of summation over spin. 
We note that actually n(k) = =n (W, (k)). The value of the function 
(11.12) enters into Eq. (6. 15’) when v=0. Obviously, x (k, 0)< 0 
if an/ow, <0. We will assume this condition is satisfied. The 
breakdown of this condition would mean that the system does 
not exist in thermodynamic equilibrium, 

According to Eq. (6. 15" ) the screened potential of an isolated 
impurity atom in the coordinate representation has the form 


tkr 
V(r) = oo \gk 2. (11.13) 
a H+ ———| Kk, 0)i 


Here p,(k) is the Fourier transform of the charge density in the 
impurity atom with its nucleus located at the origin of the co- 
ordinates. Considering the impurity atom as a point charge, we 
should set p,(k) = Ze?, 

The form of the function U(r) depends on the analytical prop- 
erties of #(k, 0). In the absence of doping the distribution func- 
tion n(k) in the present approximation is simply the Fermi function 
for renormalized energy W,(k). Thus, in the case of complete 
degeneracy (J =—0) the right side of Eq. (11.12) has a logarithmic 
Singularity at some value of k (when W, = h%k#/2m,, it is easy to 
demonstrate this by calculating the integral on the right side of 
Eq. (11.12) directly [12]; the logarithmic singularity occurs at 
k=2kr). This circumstance leads [70] to a relatively slow de- 
crease in V(r) at large separations: 


(2a, 

In the presence of an impurity the situation is altered. AS we 
saw, here there is no longer a sharply defined Fermi surface; 
consequently the function #(k, 0)is regular for all values of k, 
which is easily demonstrated by substituting into Eq. (11.12) the 
distribution function in the form of (11.2). Correspondingly, at 
large separations the screened potential of a point charge [the 
right side of Eq. (11.13) for p,(k) = Ze?] has the Debye form 


Zz. 2 
V(r) = —— exp (—r/r,), 
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where 
—~2 4ne? ,. Spe 2e? (2m,)"/# ~ 1/ on (W,) 
Ss k, ey eee i. ee ee . 
" e a eee 0) Enhs \ Uh OW, aM 
(11,14) 


The second line here is obtained by an expansion of the right 
side of Eq. (11.12) in a series in terms of k when v = 0.* 

The region where the Debye form of the screened potential 
is valid plays an important role in problems of interest to us. 
However, it is also important in calculating the parameter a 
from Eq. (8.24). By substituting Eq. (11.13) into Eq. (8.24), we 
obtain 


piss ee \4 k G oe — 5elk, o| (11.15) 


2e2n2 


The relatively small values of k’, smaller than (or comparable 
with) | (e2/ex) 4(k, 0)|, play a major role in the integral, There- 
fore the function # can be replaced by its limiting value when 
k =0, as a result of which Eq. (8.25) is obtained with the screen- 
ing radius of Eq. (11.14). Taking account of the subsequent 
terms of the expansion of # in powers of k? gives, as one can 
Show, only minor corrections. 

If the Fermi function is substituted into the right side of Eq. 
(11.14) as n(W,), then the usual expression for the Debye radius 
is obtained for a given degree of degeneracy [in particular, Eq. 
(8.1) is obtained when 7 = 0]. To take into consideration the in- 
teraction of electrons with impurity atoms, it is necessary to 
use Eq. (11.2) (for 7=0). Then the integral on the right side of 
Eq. (11.14) can be represented in the form 


—21m a —)' ( du {ip + a; (u?)} exp {ipu? + a, (u2)}, (11.16) 


where the function a, (u*)is obtained from a, (s)by replacing s by 
uw, and the prime signifies differentiation with respect to the 


*Attention should be directed to the sequence of the limiting process: first, it is 
necessary to set vy =0 and only then to make the expansion in terms of k. With an- 
other sequence for the limiting process we would obtain o(0,v) = 0. The quantity 
of (0,v) is the first term of the expansion of o% (k,v) in terms of & for a given v. A 
function of this type arises, for example, in the theory of plasma oscillations [12]. 
The difference between the two limiting expressions, lim or (k,O) and lim a (0,v), 


k — 0 vag 
is not associated with the approximations assumed above; it also remains in the exact 
theory. 
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entire argument. Here we employ dimensionless units (l= 
2iip- =), 

Since p is a large quantity, the integral of (11.16) can cer- 
tainly be calculated with the aid of approximation (9.20), Then 
we easily obtain 


wr) 
feed cos (| (1 se * (11.17) 
WL 


We note that the Debye radius (8.1) can be represented in the 
form (2/V na )(2mpuo/n2)/*, Thus the corrections to the screening 
radius because of the interaction of electrons with one another 
and with the impurity are not significant, which also confirms 
all the previous estimates. In particular, according to Eqs. 
(11.17) and (10.7) the estimate 


epee ee 
m 


used before is valid. 


SUMMARY OF THE RESULTS 
OF CHAPTER 11 


As a result of the interaction ofcharge carriers with the im- 
purity atoms their momentum distribution function n (k) differs 
from the Fermi function. In particular, neglecting quantities of 
the order of m/M, where M is the mass of an impurity atom, an 
expression is derived for n(k) that is continuous, together with 
all its derivatives for any temperature, including 7=0 (the Fermi 
‘‘step”?’ is *‘smeared out’’), Consequently, the concept of the 
Fermi surface in a heavily doped semiconductor has no unam- 
biguous meaning. This leads to a relatively simple form of the 
screened potential produced by the individual impurity atoms. 
Thus, at asymptotically large distances from the impurity atom 
the standard Debye equation is valid. The screening radius is 
only insignificantly different [see Eq. (11.17)] from the Debye 
value for a degenerate gas. 


CHAPTER 12 


Static Electrical Conductivity 


From an experimental point of view a study of the kinetic 
coefficients which characterize the behavior of a material in d-c 
and varying fields—the complex electrical conductivity, the 
Hall constant, thermoelectric power, etc.—is of the greatest 
interest. 

In the present chapter we will consider the problem of the 
static electrical conductivity of a heavily doped semiconductor. 
We restrict ourselves to the case of quite low temperatures 
when scattering by a charged impurity, but not by phonons, 
plays a role. One can imagine that the concept of ‘‘quite low’’ 
temperature encompasses quite a broad range of 7 in a given 
case—in Ge, for example, up to and including room tempera- 
ture. It is known from experiment (see, for example, 
[71, 72]) that the electrical conductivity of n- and p-type Ge 
containing 101% -10'° cm-3 atoms of the doping impurity is ap- 
preciably smaller by a factor of 3 to 10 than the electrical con- 
ductivity of lightly doped specimens even at room temperature. 

The characteristics of a heavily doped semiconductor, con- 
sidered in Chaps. 2 and 3, make it somewhat hazardous to cal- 
culate its electrical conductivity with the aid of the standard 
kinetic equation method and one is forced to turn to the general 
equation (5.23). The problem thus involves a calculation of the 
two-particle Green’s function, Eq. (5.24). 

Later on we will make use of the ‘‘one-band’’ effective mass 
approximation with the simplest quadratic dispersion relation. 
Then Eq. (5.23) assumes the form 


89 
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Tap (X) = — Slim = (ck dks dks (Rrekan + Raabe): 
se-2nt(s, keke) 36 (kx, key ks kus v), Can 
where 
Env 
and 
6 (ki, Ke, kg, ka; v) = \ dx dy dx’ dy’ exp 2mi ((kix— key-+-k,x’ — 
— kay’} H(x, y; x’, y’; v). (12.2) 


Here the Fourier transform is generated independently in all 
coordinates because before the averaging over the impurity con- 
figurations the system is not yet spatially homogeneous. After 
the averaging we obtain, obviously, 


H (ki, ke, kg, ka; v) = 8(ky + ky — ky — ky) 4" (ki, ke, kg; v)- 
(12.3) 


The conductivity tensor (12.1) no longer depends on the coordi- 
nates. Introducing the symbol 


x’ (k, k, k's v) = R(k, k's v), (12.4) 
we Obtain from Eqs. (12.1) and (12.3)* 


ie tim m7 \ dkd k’kakpR (k, k’; v). (12.5) 


m? Eo 


02.8 = — 


As we have seen, under conditions of extremely heavy dop- 
ing (A<!), taking the Coulomb interaction into account simply 
reduces to using the screened potential and toa renormalization 
of the dispersion relation [replacement of W(k) by W,(k)]. Before 
averaging over the impurity configurations we have a problem 
with the formally additive Hamiltonian, and the function #(x, y; 
x’, y’; v) is given by Eq. (11.9). The auxiliary function g appear- 
ing here, which is determined from Eq, (11.10), is easily ex- 
pressed in terms of the exact one-particle Green’s functions G, 
and Ga. We set (before averaging !) 


*Here, and later, we are following [73]. A method similar in concept was devel- 
oped in [74]: however, the approximation used there apparently does not provide a 
correct limiting process for the case of the kinetic equation. 
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Gra (x, > ae v) =; \dkd k’erni(ha—k'x’) G,,a (k, k’; v). (12.6) 


Taking Eqs. (5.9) and (5.11) into consideration, weeasily derive 
(12,6) from Eqs. (11.10) and (12.2). Now 


+00 
# (ki, ky, ks, ky; v) = \ dv’ [1 --e-Bihy’—w))-1 (G(ky, ky; v’—v)- 
-[G, (Ke, k3; v’) — Ga(Ke, kg; v’)] + G, (ke, ks; v’ + v)- 
-[G, (Ky, Ki; v’) — Ga(Ka, ki; v’))}. 
(12.7) 


Substituting Eq. (12.7) into Eq. (12.1) and integrating by 
parts, we can obviously perform the limiting process F-0. 
Averaging the result over the impurity configurations and taking 
account of the identity [12] 


G,(k, k’; v) = Ga(k, k’; v), 


we obtain 
TO ape 
lim °- R(k, k's; v) = —i ‘Cage 24 ne Kk, ky v)— 
E+») OE 
— Reke(k, k’; vj}, (12.8) 


where the functions A; and K, are defined by the equations: 


ee ee 
Gk, k's \) Ga lk", ky v= 8K — Rk) K(k, ky vy. 49) 


—_—— 
G, (k, k’; v)G, (k", k; v) = 8(k’ —k”) Ko (k, Kk’; v), (12.10) 
and 
np (E’) = [e2’—) + 1. (12,11) 


The summation over spins again gives a factor of 2 and we ob- 
tain finally 


0 
ee (ae et dk dk’, ke {Ki(k, k’3 Vv’) — 
m2 


— Re Ky (k, k’; v’). (12.12) 


We note that in a cubiccrystal, where the electrical conductivity 
tensor degenerates to a scalar, only the part of R which is 
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symmetrical with respect to a transposition of k and k’, is of 
interest. Moreover, in an isotropic system the integral with 
respect to k’ in Eq. (12.12) is obviously proportional to kg: 


\ dk’kp { ... =F (Rk, Vv). (12.13) 


On the other hand, the usual equation for the electrical conduc- 
tivity that results from the kinetic equation under the conditions 
being considered has the form 


26th? =e —, 


| dk ka ket (Ws) (12.14) 


Oap == 


where t is the relaxation transport time. Substituting Eq. (12.13) 
into Eq. (12.12) and comparing the result with Eq. (12.14), we 
See that the expression 


ae 
\ dE’ - 


plays the role oft dn,w/dw, In particular, for T =0 we have 


— 


F (k, E’) (12.15) 


F i(k, p) = t6(W, — 1). (12.16) 


From this it is clear that the usual concept of relaxation time, 
as a regular function of W,, can be introduced only if the func- 
tion F(k, jt) is a delta function. This is precisely the case in the 
first nonvanishing approximation in terms of the small param- 
eter A. 

According to Eqs. (12.12), (12.9) and (12.10), the static 
electrical conductivity problem reduces to a calculation of the 
average product of the two one-particle Green’s functions. This 
quantity, of course, is by no means the same as the product of 
the average functions calculated, for example, in Chap. 7 or 
Chap. 8. Therefore it is necessary to perform the calculation 
anew, to a certain extent. The results obtained earlier do, 
however, give some orientation. According to Chaps. 9 and 10, 
when 2<1 the electron density in the ‘‘tail’’is relatively small, 
as is also the correction to the Fermi level. For this reason 
one can look forward to obtaining reasonable results simply with 
the aid of perturbation theory of the type developed in Chap. 7. It 
Should only be kept in mind that the principal term in the elec- 
trical conductivity will obviously be of the order of1/n. There- 
fore, a direct expansion ofo in a power Series in terms of the 
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Fig. 6. Equation (7.187), 


density is meaningless. The perturbation theory, however, can 
be used to formulate the integral equation which the functions K, 
and K, (see [46]) should satisfy. This is precisely the prob- 
lem that was presented in [54] in which the impurity re- 
sistance of metals was calculated by the methods of field theory. 
In essence, in the present case it is only necessary to take ac- 
count of the corrections discarded in [54] associated with 
the effect of the impurity on the energy spectrum of the charge 
carriers and to allow for multiple scattering. This program 
has also been carried out in [73]. Just as in [54] it is 
convenient to make use of a graphical technique by comparing 
the unperturbed Green’s function (dashed line), the perturbed 
function (solid line) and screened potential (wavy line). Then 
Eq. (7.18) will be represented by Fig. 6, and the second-order 
terms in Eq. (7.20) by Fig. 7 (these are precisely the terms 


Ce the 1s a, es et 
Fig. 7. Second-order terms in Eq. (7.20). 


that were explicitly calculated in Chap. 7). Multiplying the two 
series of Fig. 6 and averaging the product, it is easy to verify 
that to a first approximation the function ,; is defined by the 
equation of Fig. 8 (the same is also obtained for K,). The solu- 
tion of this equation is of order 1/n. The corrections are derived 


Fs =F 


Fig. 8. Integral equation for the function K, (first approximation). 


from the intersections of the wavy lines in a diagram such as 
Fig. 9. We note that lines appear in Figs. 8 and 9 which corre- 
spond only to the exact (‘‘dressed’’) one-particle Green’s func- 
tions. In calculating the corrections from perturbation theory 
this means that the mass operator M+, Eq. (7.21), should be 
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calculated taking account of not only the diagrams of Fig. 7 but 
also Figs, 10 and 1l. 


SESS eR, PS ,.. 


Fig. 10. Correction to the mass operator M,. 


sl ce thes 


Fig. 11. Corrections to the 
mass operator \f,. 


For the explicit calculation in [73] the Debye potential (7.35) 
with the screening radius (8.1) was used. (We have seen in 
Chap. 11 that this approximation is valid). For the electrical 
conductivity at 7=0 there is derived 


eae (ae ee) EDU + BiE)}, (12,17) 
322m |) dE 7 


where t is the usual transport relaxation time; the function k(E) 
is the solution of the equation 


E—W,(k)=0, 


and the function B(F) has the form 


BR Sy (ee eZ 
(R) 13a5, 73 (442 + 72) fer 7 ° | 12.18) 
(12. 


Here V and Z are slowly varying functions of kro, Their explicit 
forms are quite unwieldly; a numerical estimate is introduced 
below. 

Obviously, the first term in the braces in Eq. (12.17) gives 
the usual result obtained with the aid of the kinetic equation 
(see, for example, {37]); the function B(E) describes the 
corrections associated with two-stage scattering and the effect 
of the impurity on the charge-carrier energy spectrum. 
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When 7=0, the integration in Eq. (12.17) is performed di- 
rectly and we obtain 


o =o, — Ao, (12.19) 


where oo is the usual expression [37] 


6, = elke). (12.20) 


m 


The correction As/s, is negative and oforderA. With parameters 
characteristic of germanium, it amounts to about 25%. A similar 
result is also obtained at a finite temperature. 

It should be noted especially that the quantity Ac/s, exhibits 
an extremely weak dependence on the impurity concentration 
under the assumption of complete degeneracy of the electron 
gas. For this reason the results discussed here have not yet 
been checked by experimental data [71, 72]. One might have 
assumed that the difference between the actual screened potential 
and the Debye potential plays a definite role. Actually, the 
Debye potential corresponds to the asymptotic form of Eq. (6.15) 
for k = 0. A calculation with the total polarization operator [73] 
showed, however, that the corresponding correction amounts to 
about a tenth of a percent. Thus the question of a valid interpre- 
tation of the results of [71, 72] remains open.* In connec- 
tion with this attention should be directed to the fact that for the 
concentrations stated in [71, 72] and the parameters char- 
acteristic of germanium, the asymptotic form of ~< 1! is not 
attained. Under such conditions the equations derived above 
might be altered for a purely impurity scattering and, moreover, 
another relaxation mechanism—for example, the scattering of 
electrons taking account of transfer processes—might be im- 
portant. These questions have not yet been properly examined. 


SUMMARY OF THE RESULTS 
OF CHAPTER 12 


The expression for the static electrical conductivity of a 
heavily doped semiconductor, associated with scattering at a 
charged impurity, does not differ greatly from the standard 


*Footnote added during proof. During the printing of this book it was learned that 
the results of [71,72] are in good agreement with normal kinetic theory taking account 
of the anisotropy of the equi-energy surfaces [119]. 
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result of kinetic theory. In the final analysis this is explained by 
the small effect of the impurity on the charge-carrier energy 
spectrum in the region of energies close to the Fermi level. 
The correction to the result of the kinetic equation, calculated 
for the case of complete degeneracy, exhibits an extremely 
weak dependence on the impurity concentration, 


CHAPTER 13 


Absorption of Electromagnetic Waves 
in a Heavily Doped Homopolar 


Semiconductor 


The optical properties of heavily doped semiconductors have 
in recent times been the subject of rather extensive experi- 
mental research. 

Let us first consider the absorption associated with band-to- 
band transitions. 

In relatively recent papers [75-80], referring to the inter- 
metallic compounds n-InSb, n-InAs and p-GaAs, it has been 
shown that the threshold for band-to-band transitions is shifted 
toward short wavelengths. This effect was interpreted [78] as 
the result of the degeneracy of the electron gas: in neglecting 
the effect of the impurity on the band structure and for p>0, 
tssT, the threshold frequency thr is defined, obviously, by 


othr = A+ LL (13.1) 


(the zero of energy, as always, coincides with the bottom of the 
conduction band). 

A more detailed study of the absorption edge was carried out 
in later papers [81-86]. Experimental data for Ge [82-86] can 
be more or less satisfactorily interpreted with the aid ofa 
hypothesis concerning a narrowing of the forbidden band with 
an increase in the impurity concentration (both donors and ac- 
ceptors). Numerical results for the amount of narrowing in 
different papers are somewhat different—from 0.03 eV to0.09 eV 
(in n-Ge with ny ~4 x 10'9 cm-3). In Si, according to [81], 
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the effect does not exceed 0.03 eV. Itshould be noted that a study 
of the shape of the ‘‘band-to-band’’ absorption bandin the region 
close to the threshold meets with many experimental difficulties, 
That is, absorption in this frequency region is associated not 
only with band-to-band but also with “‘intraband’’ transitions 
(absorption by free carriers). In view of the high concentration 
of charge carriers the effect is by no means small. In order to 
avoid this, the absorption is studied in the frequency region hbe- 
low the threshold when band-to-band transitions are impossible. 
Having determined in this manner the frequency dependence of 
the coefficient of light absorption by the free carriers, itis 
extrapolated to the region of higher frequencies. Subtracting the 
result from the experimental value of the absorption coefficient, 
values are obtained which describe the band-to-band transitions. 
The latter, unfortunately, is almost an order of magnitude 
smaller than both the minuend and the subtrahend. Therefore 
the above procedure cannot claim a high accuracy. Moreover, 
in the region of energies close to the Fermi level (in an uncom- 
pensated specimen) the ‘‘bare’’ dispersion relation, determined 
experimentally for states near the bottom of the band [114], can 
no longer be exact. Consequently, not too much importance 
should be attached to the divergence inthe results at the present 
time. However, an effective narrowing of the forbidden band 
does, in fact, obviously occur in heavily doped specimens and 
one should attempt to interpret it from the point of view of the 
representation of the band structure of such systems previously 
discussed. 

To calculate the absorption coefficient a, associated with 
band-to-band transitions, let us use the familiar relationship 
between a and the real part of the electrical conductivity ata 
given (angular) frequency wo: 


dics = Reg (a) (13.2) 
cVe 


Here c is the velocity of light in vacuum and e is the real part 
of the dielectric constant. In the absence of localized levels the 
value of e in a semiconductor is determined primarily by the 
crystal lattice; in the frequency region of interest to us the 
value of e() (in a homopolar crystal) is practically constant. 

For cubic crystals and isotropic media Eq. (13.2) is derived 
directly from Maxwell’s equation under the condition that the 
frequency » exceeds the plasma frequency of the electron gas 
and e > (4:/S) Reo(w), In the infrared region these conditions 
are usually satisfied. It is significant that Eq. (13.2), deriving 
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from a phenomenological theory, is in no way associated with 
any model representations (e.g., with the hypothesis concerning 
the existence of the dispersion relation). This makes it possible 
to use it in the problem of interest to us when the standard 
method of the theory of quantum transitions is not sufficiently 
conclusive in view of the characteristics of a heavily doped 
semiconductor. 

To calculate e(w) we make use of Eq. (5.22).* Let us turn to 
the concept of second quantization, having selected 


a(x, S) = So(s) exp (— px Ue (X) (13.3) 


as the basic system of Bloch functions describing one electron 
in an appropriate ideal lattice. Here S is the spin function, s 
and o are the coordinate and spin quantum number; the index A 
designates the set of quantities {p, /, o}. 

The current density and velocity operators will have the form 


= 5) Aljald’ Jara, (13.4) 
rd 


A, 
= (Aug | dana, 
AA’ 


where (A\ja'A’) and (A|ug|A‘) are the matrix elements of ja and vu, 
in the system of functions (13.3). In the absence of a magnetic 
field and spin-orbit coupling they have the form 


(RL ja| R= Bare (P) + (1 bir) B0Spp0ir'(p), (13.5) 


(4 |p 1%’) =633.0p (P) ++ (1 811) Ba0-Spp-tiv- (P), (13.6) 
where 2, as before, is the volume of the system, 
v(p) = VW 7(p), v7 (p) = — = =| dx re lr - (13.7) 
) 


and i) is the true electron mass. 
By definition, when /’ J, 


¥ 
> = | Ore" (p) |? wae fae (13.8) 


2M 


*We follow [87] here. 
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where /fi:(P) are dimensionless parameters—the oscillator 
strengths for the /’=/ transition. These quantities are the 
essence of the spectral characteristics of an ideal crystal and 
in the problem under consideration they should be assumed to 
be known. Actually, they can be reduced from optical data only 
with difficulty (and somewhat ambiguously). For now, however, 
it is sufficient to know that they are smaller than unity by an 
order of magnitude and do not change rapidly with a change in p. 
The first of these statements is a consequence ofa familiar 
sum rule (see, for example, [88]); the Second follows simply 
from the fact that f is dimensionless. Being a dimensionless 
scalar quantity characterizing an ideal lattice, the oscillator 
strength can depend only on the combination pd/A, where d is 
the lattice vector. 

According to Eqs. (5.22) and (13.4) the complex electrical 
conductivity is expressed in terms of a two-particle Green’s 
function of the form 


KH ,,,| 0,0, > = 0()0 <la, (Nar, (8), &, (0) a, (>. (13.9) 


In general, to calculate this it is necessary to solve an inte- 
gral equation (as was done, for example, inChap. 12 for the case 
« =0). In the present section, however, we are only interested 
in a sufficiently high frequency, iw~A, in which the indices A, 
and 42 (as wellas 43 and 24) Should refer to different bands (band- 
to-band transitions). It can be demonstrated* [87] that in this 
case the approximation 


<KA;,0}, | An,00, > 2 ¢ da, (t).2,(0)> G (Aa, Ages B-E 


-- €ay, (QO) aa, (f)> Guha, Ari 2) 

is valid asymptotically for large values of \. 

The quantity !/(A\1)*, where t is ofthe order of the relaxation 
transport time, plays the role of the small parameter that en- 
sures (in zeroth approximation) the validity of Eq. (13.9’). The 
approximation (13.9 ) cannot be used in the consideration of in- 
traband transitions when the finiteness of t is of prime im- 
portance; for band-to-band transitions, however, this approxi- 
mation is sufficient. 

Taking into consideration the relationship (5.9) [and Eq. 
(5.11)], it is easy to obtain from Eq. (13.9) 


(13,9’) 


*See also [89, 90] where an equation similar to (13.9) is derived with the aid of 
causal Green’s functions. 
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Im < ay,03,| 09,0), Sy 2 \ dv’IMG, (Az, X33 W’). 
{Im G, (Ay Aq; v’ — v)ae(v’— v)— [inG, (Ayay3; vv! — v)ne(v')}, 
(13.10) 


where the symbol «...», represents the Fourier transform of the 
function (13.9) and nr is the Fermi function 


np (v) = [ebay Yr, (13.11) 


In carrying out the averaging of Eq. (13.10) over all impurity 
configurations, it should be noted that the asymptotic approxi- 
mation (13.9°) is also valid for the averaged Green’s functions. 
The latter, when ‘interaction between bands’’ is neglected, 
have the form 


Gh, M5 v) = OarG(As v). (13,12) 
Taking this relationship into consideration, substituting Eq. 


(13.4) into the general equation (5.22) and making use of the 
approximation (13.10), we obtain 


Sire? e 
Re dap (#) = —— : Dvir) viu(P) » 


oe 


x \ dv’ Im G, (p, é; v’)ImG,(p, d; v’ — v){n (v’ — v) — ne (v’)} 
(13.13) 


for the real part of o.g associated with band-to-band transitions. 

In a cubic crystal this tensor degenerates to a scalar; accord- 
1 ef y) 

ingly the product u7-u;, can be replaced here by > >) vi [*bzs. 


Taking into account the definition of Eq. (13.8) and changing from 
the summation over p to an integration, we finally obtain 


e2 
SS SSS d , Ay y 
a(0) = ae \ae D) feu(P) Ars (p) 
+00 
x \ dv’ ImG,(p, l’, v’) ImG,(p, lf; v’ — v) {ne (v’ — v) — ne (v’)}. 
oo (13,14) 


Here wo =2nv has been used and the quantity 
Wy (p) — Wi (p) = Arve (P) (13,15) 
has been introduced. 


102 ELECTRON THEORY OF SEMICONDUCTORS 


Obviously « is an even function ofw. Later we will assume 
w>O. 

Clearly, A; ,(p) is the separation between dispersion curves 
(in an intrinsic material) at a given point of the Brillouin zone. 
The fact that precisely this quantity enters into Eq. (13.14) does 
not mean, however, that we are considering only direct transi- 
tions. Actually, the Green’s function entering into Eq. (13.14) 
is exact. In principle it also takes account of the interactions 
responsible for indirect transitions. This appears formally in 
the fact that its imaginary part is not generally a delta function. 
On the other hand, the energy condition 


Ayer (p) = Ay, 


which distinguishes just direct transitions, is obtained only if 
one sets 


ImG,(p, l: v) = 6 (v - ——W, () | (13.16) 


By substituting Eq. (13.16) into (13.14) we obtain the standard 
equation for the absorption coefficient in the case of direct 
transitions: 


X (0) —— re \ dP S) fire (9) 8 (hv — er (0) » 


2mych? Ee 


W 
of rr ea) —N"rF We /h)} (13.14 ) 


When the interaction of electrons with impurity atoms, 
phonons, etc., is taken into account and when Eq. (13.16) is no 
longer valid, the terms in the sum over /, /'in Eq. (13.14) can 
be divided into two groups. In the first of these the argument wv’ 
corresponds to the higher band (W,. >W,); in the second group 
the same is true of the argument v’—v. When the concept of 
different bands retains its meaning in the presence of an im- 
purity, the terms of the second group can be neglected: the 
‘‘smearing out?’ of the band levels is not sufficient that the 
states deriving from the ‘‘former’’ /th and /’th bands exchange 
places. Only the case of transitions between different ‘‘sub- 
bands,’’? comprising one degenerate or nearly degenerate band 
(for example, the valence band in germanium), can constitute 
an exception. In such systems the concept of band-to-band 
transitions in the presence of an impurity generally loses a 
definite meaning, 
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Taking into consideration what has been said, let us now 
specialize Eq. (13.14) for the case of two simple bands—the 
valence (v) and the conduction (c) bands. In this situation the in- 
dices /’ and | assume the values c and v, respectively, and we 
obtain 


e2 
es ea er \ dP feo (p) Aco (P) X 
+00 
x ( dv’ ImG, (p, c; v’) ImG, (p, v; v’ — v) {n- (v’ —v) —ne (v’)}. 
ase (13.17) 


Equation (13.17) is derived under the assumption that the 
system of electrons exists in a condition of static equilibrium 
in the absence of a radiation field. One special nonequilibrium 
case, namely, when energy equilibrium and crystal momentum 
equilibrium in each band has already been established but the 
densities of electrons and holes are not equilibrium values (in- 
jection), can be of experimental interest. It is well known that 
such a situation is relatively easy torealize under nonstationary 
or stationary conditions. The reason for this is the slowness of 
the charge-carrier recombination process compared with other 
relaxation processes. 

Under injection conditions we have not one but two independ- 
ent conserved quantities—the conduction electron and hole den- 
sities. Consequently, in describing the system it is necessary 
to introduce two quasi-Fermi levels, jte and w»[91}. It is easy 
to show that in this case Eq. (13.17) must be rewritten in the 
form 


+o 
ee —\dpfeo(p) Ac (p) \ dv’ ImG,(p, c; v’) 4 


mm cwht Ye « 
—0Oo 


hv’—hvtAtu =r 
x ImG, (p, 9; vs {l exp" 7 maa] = 


Av’ — |, —1 
—[exp—* 41] 3 
(13.17°) 


Here p, is measured from the bottom of the conduction band and 
uw, from the bottom of the hole band (i.e., downward from the 
top of the valence band). In particular, when 


He +A + pa > Av (13,18) 
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the absorption coefficient becomes negative. This obviously 
means that the system does not absorb but emits radiation. The 
emission coefficient is given by the absolute magnitude of a; the 
inequality (13.18) is nothing more than the ‘‘inversion condition 
at frequency v,’? expressed in exact terms without using any 
‘fone-electron’? concepts. 

As seen from Eqs. (13.17) and (5.15’), the absorption coeffi- 
cient is not generally expressed directly in terms of the density 
of states in the upper and lower bands. A measurement of a at 
an arbitrary temperature can only give rather indirect informa- 
tion which is difficult to interpret concerning the energy spec- 
trum of a heavily doped semiconductor. The situation is altered, 
however, at sufficiently low temperatures that satisfy the condi- 
tion T<pn (to be specific, we are considering an n-type speci- 
men with the Fermi level, as always, measured from the 
bottom of the conduction band in an intrinsic material). In this 
case nr (v’)= 6(u —Av’) and the integral with respect to v’ in Eq. 
(13.17) is taken over the limits from p/h to n/h+v. AS we saw 
in Chap. 9, near the Fermi level the impurity has a relatively 
weak effect on the energy spectrum and the density of charge- 
carrier states. The same also applies to the effect of the im- 
purity on the ‘‘foreign’’ band. In neglecting these small effects 
here we formally go over to the usual theory of light absorption 
due to band-to-band transitions.* The only difference lies in the 
fact that the forbidden band gap is now determined by the ener- 
gies, ‘‘renormalized’? due to the Coulomb interaction, of the 
electron and hole: 


A—> A, = Min {(W,.(p) — W, ov (p)}. (13.19) 


By altering the dispersion relation, the Coulomb interaction 
can in principle also shift the extremum from that point in the 
Brillouin zone which it attains in an unperturbed system. An 
analysis shows, however, that in the case of spherical or 
elliptical energy surfaces this does not occur. Thus, according 
to Eqs. (13.19) and (7.12), 


Aa (13,19’) 
Where “.and VM; are the mass operators for electrons and holes, 


taken at corresponding extremal points (for example, forp =0). 
In a heavily doped n-type specimen the quantity M,, turns out to 


*An indispensable summary of the results of this theory can be found, for example, 
in {g6] and [92]. See also [93]. 
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be very small, as is easily shown. Consequently, for the sim- 
plified dispersion relation (7.13) and T = 0 we obtain from Eqs. 
(13.19’) and (7.14) (neglecting the ratio x/k-) 


he eee —2{ SP eer x oe (13,20) 
It E E 

In the conduction bands of germanium and silicon the equi- 
energy surfaces near the bottom of the bandare a set of several 
ellipsoids, The centers of the latter (minimum energies) lie on 
the symmetry axes of the reciprocal lattice at equivalent points 
of the Brillouin zone. The dispersion relation near each min- 
imum has the form 


em 
W (p) = = 


+ se » My <M. (13.21) 


Here the 0z axis is directed _ the cies axis of the ellipsoid 
and the components of the crystal momentum are measured 
from the minimum point. In this case a simple calculation from 
Eq. (7.7) gives 


Piha 18 (2)" UE 5 Me eee, (iG 22) 
nt a Vxty) Vy 
where ¥ = m,/m, and z is the number of minima. 

Equation (13.22) agrees qualitatively with the results of 
[82-86]. Numerically, for germanium we have [94] « = 16, m, = 
1.6 my, mt = 0.082 m,, 2= 4, For na = 4X 10! cm® this gives, 
according to Eq, (13.22), 


A, —A--—Q0.1 eV. 


This quantity should be compared with the data mentioned at the 
beginning of this section. A discrepancy by a factor of 1-1/2-2 
Should not be taken too seriously, both in view of the poor ac- 
curacy of the experiment and also becausSe the ‘‘small parameter 
47? is close to unity, and Eq. (13.21) is only suitable for order- 
of-magnitude estimates. 

When the effect of the impurity on the density of states is 
taken into account, the result is complicated and generally no 
longer reduces simply to a narrowing of the band in the literal 
sense of the word. The frequency dependence of « is also altered; 
however, as already stated above this effect isrelatively small. 
If, nevertheless, it is taken into account, then the concept of 
band narrowing loses definite meaning: the ‘‘shift??A,—A in this 
case is some effective value that approximately describes the 
actually much more complicated spectrum variation. With suit- 
able accuracy, experiments of different types carried outon the 
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Same specimen should generally give different values for the 
effective band narrowing. The absence of such discrepancies 
would most likely indicate simply the loss of all ‘‘impurity”’ 
effects in view of insufficient experimental precision. 

To reduce the role of absorption by free carriers, experi- 
ments with heavily compensated specimens are essential. A 
study of such materials is also of interest in the sense that with 
sufficiently precise compensation the Fermi level can fall to 
the middle of the forbidden band. In this case, according to Eq. 
(13.17), a direct study of the ‘‘tails’’ ofthe density of states be- 
comes possible. Such experiments are described in [108, 
110]. These results apparently agree with the concept of a “‘tail;’’ 
however, a quantitative interpretation of these data is still 
lacking.* 

Experiments on recombination radiation (under conditions of 
injection) are of vital interest in studying the effect of the im- 
purity on the spectrum of the charge carriers. When p,.>>7, and 
U,s>T, we have, according to Eq. (13.17), 


MW p/h 
i 
}a1~—\ dpfeo(p)da(p) | dv’ ImG,(p,  v)x 
Atin 
— +y 
xImG,(p, 0; v’ — v). (13.23) 


Clearly the appearance in the emission spectrum of frequencies 
less than 4,/h would be an exact proof of the existence of a ‘‘tail’’ 
in the density of states. Of more interest, however, is the possi- 
bility of investigating the form, itself, of the density of states 
in the ‘‘tail?? with the aid of recombination experiments. 
Strictly speaking, for this the Green function appearing in Eq. 
(13.23) should be calculated with the aid of one model or another 
(for example, the approximation in Chap. 8). An estimate of the 
results can, however, be obtained by noting that in the “‘tail’’ 
region the imaginary part of the Green function changes very 
rapidly; therefore the product fo.(p)4.(p) can be taken outside 
the integral sign as some constant. Moreover, the effect of the 
impurity on the ‘‘foreign’’ band is relatively insignificant and, 
consequently, by considering an n-type specimen Eq, (13.16) 
can approximate the imaginary part of the hole Green function, 
Thus, for hv<A 


*One effect - the change’in the nature of the indirect transitions associated with 
the increased role of impurity scattering compared with phonon scattering - is estab- 
lished in [108]. 
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{ 1 A 
ja| ~ —{ dp Im G, (, Cc; Vv — 7 W,o(p) — 7) Gn — W,» (p)). 
(13,23) 


Here W,,(p) is the renormalized hole energy measured from its 
minimum (in the isotropic approximation W,, = p’/2m,,). More- 
over, let us restrict ourselves to materials like GaAs in which 
the extrema of the conduction and hole bands lie at one and the 
same point of the Brillouin zone. Then the function 6 under the 
integral sign in Eq. (13.23') can be replaced by unity when 


Un > E,, (13.24) 


where E, isaconstant, introducted in Chap. 9 and defining the ef- 
fective depth of the ‘‘tail.’? Actually, under the conditions stated 
the integral expression in Eq. (13.23’) is already essentially 
terminated when W,, is of the order of a few E,. Moreover, let 
us assume that the effective mass ofa hole, m,,;, is appreciably 
greater than the effective electron mass, m,;: 


Mp > Mr. 


Then the variation of the function G, (p,c; v —W,)(p)/h —A/h) with 
the crystal momentum will be determined, primarily, by its 
argument; however, in the integration over p it can be assumed 
constant. Then, according to Eqs. (13.23’) and (5.15), 


|x| ~— pe(hv — v), (13.25) 


where 0 is a constant. of the order of A. 

It should be remembered, however, that Eq. (13.25) is only 
an estimate; coefficients of the type v7, where g is a small 
number, can be ‘‘smeared’’ in the assumed approximation. 

A detailed investigation of the absorption associated with di- 
rect band-to-band transitions in semiconductors like germanium 
could have special interest for our problem, In an uncompen- 
sated heavily doped specimen the Fermi level lies above the 
bottom of the conduction band but certainly by no means above 
the secondary minimum located in the center of the Brillouin 
zone. Consequently, one could obviously study the effect of the 
impurity on states with small crystal moments. Experiments 
of this type are described in [84, 95, 111]. To a consid- 
erable extent their interpretation is free of the difficulties men- 
tioned above because in the ultraviolet region, with which we 
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are dealing here,* the role of the free carriers is not so im- 
portant. However, the range of concentrations encompassed in 
the specimens cited was apparently sufficiently broad that the 
effect of the passage of the Fermi level through the correspond- 
ing secondary minimum was detected. On the basis of the ex- 
periments of [84, 95, 111] it can only be said that at points 
of p-space relatively remote from the extrema of the conduction 
and valence bands, the relative location of these bands is only 
Slightly altered during doping. Thus, in p- and n-5Si (doping 
elements—boron, phosphorus and arsenic) the relative shift of 
the bands does not exceed 0.02 eV [111]. No adequately detailed 
theoretical study of this problem has yet been made. 

Let us now return to the optical characteristics associated 
with intraband transitions. The absorption by free charge carri- 
ers was studied experimentally in [96]. The authors of this 
paper concluded that the dependence of the absorption coefficient 
on the wavelength 4 is parabolic (a~A?); however, in view of the 
smallness of the range of variation of 4, this conclusion is ob- 
viously not definite. 

An approximate theory of this type of absorption was given 
in [115, 116]. The authors restricted themselves to fre- 
quencies appreciably greater than the collision frequency. In 
this case the problem essentially reduces to a study of the ab- 
sorption of high-frequency electromagnetic radiation in a de- 
generate plasma (see, e.g., [117]); the most subtle character- 
istics of a heavily doped semiconductor are not mentioned here. 
A satisfactory theory for the case of longer wavelengths ap- 
parently does not exist at the present time (mid-1963). 

The contribution of free carriers to the real part of the di- 
electric constant, Ae, has been measured (-Ge in [96], n- and 
p-Ge and Si in [97]). Measurements ofthis type make it possible 
to determine the ‘‘optical’’ effective mass m* according to the 
familiar relationship 


4nn,e* 


Ae — (13.26) 


- 
Equation (13.26) is valid for wtS +1, where t is the relaxation 


transport time. The validity of this inequality can be established 
experimentally simply from the conclusion that the curve Ag (w?) 


*In actuality the experiment involves a measurement of the frequency dependence 
of a, the reflection coefficient in the ultraviolet region. As is known, the absorption 
coefficient, in which we are interested, can be directly deduced from the reflection 
coefficient when the dispersion curves are known. 
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is a parabola. According to [96] the quantity mis slightly 
dependent on the impurity concentration up to n~ 5 x 10” cm, 
This result is not especially astonishing, since, as one can show 
[12, 98], the optical effective mass, defined by Eq. (13.26), is 
exactly the same as the bare mass when the unperturbed dis- 
persion relation W(p) remains quadratic in the energy region 
being studied. The relationship 


Ree (o) = — “Ima(o) 


follows directly from Maxwell’s equation, so, according to Eq. 
(5.22), 


Im (at < ja(t), 0g (0) > et. (13.27) 


It is apparent from physical considerations that with the onset 
of the relaxation mechanism the correlation function < [/a(é), 
ve (0)]. > decays with time. It follows that the integral in Eq. 
(13.27) decreases with an increase in w, and the contribution in 
the parabolic function Ae(w—*) gives only the firstterm: m* =m, 
Moreover, comparing Eqs, (5,21) and (5.22), we obtain 


: 9 
ift .€ 


IL wii 
bap cama 77 Abre ds]_\, (13.28) 


mG 


where the operators j. and d,; are taken at the same instant of 
time. For a quadratic and isotropic dispersion relation,* W (p) = 
p?/2m , the commutator in Eq. (13.28) is easily calculated: 


° ° ne” 
([ja, ds], = — th Oa, 


and, consequently, m=, independent of the nature and mag- 
nitude of the potential energy of the electrons. In the case of 


ellipsoidal equi-energy surfaces, when W (p) = > p, 2m we 
obtain in a Similar manner as 
1 l { 
gg ot ee 13.29 
m 3 2 my ( 


*Let us recall that Eq. (5.22) has meaning only when used within the framework of 
the effective mass method: the unperturbed system is assumed to be spatially homo- 
geneous (in particular, the electron density 1s assumed to be constant). 
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Corrections associated with the interaction of the charge carriers 
with the impurity atoms, with each other, etc., may be made to 
this. 

For a nonquadratic dispersion relation the situation is al- 
tered. Thus, if the group velocity is given by the expression* 


V (P) = VW (p) = pf (0’), (13.30) 


where f is an arbitrary function, then a calculation of the com- 
mutator in Eq. (13.28) yields [12]: 


<= (apni lie) +— pt wr}. (13.31) 


Here the prime designates a derivative with respectto the com- 
plete argument; n(p) =n(—p)is the momentum distribution func- 
tion, normalized to unity: 


\dpn(p) = |, (13.32) 


In the special case when f = 1/m = const, Eq. (13.31), together 
with Eq. (13.32), yields, of course, the previous result. If, 
however, the function [ does not reduce to a constant, then the 
integral in Eq. (13.31) can no longer be calculated simply with 
the aid of the normalization condition and the result will depend 
both on the specific form of f(p?)and on the form of the distribu- 
tion function, i.e., ultimately, on the impurity concentration, 
the temperature, etc. 


SUMMARY OF THE RESULTS 
OF CHAPTER 13 


The narrowing of the forbidden band ina heavily doped semi- 
conductor, having been observed experimentally in a number of 
cases, can be explained by the Coulomb interaction between 
charge carriers [Eq. (13.19), (13.21)]. The effect of the inter- 
action of electrons with the impurity on the band-to-band absorp- 
tion spectrum cannot be described only with the aid of a narrowing 
of the band—the form of the frequency dependence of the absorp- 
tion coefficient itself is altered. However, in uncompensated 
specimens this effect is relatively small. 


*In the case of n-Ge type anisotropy, Eq. (13.30) must be written separately for 
each ellipsoid. 
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A measurement of the frequency dependence of the recom- 
bination radiation coefficient can in some cases give direct 
information concerning the shape of the density of states in the 
forbidden band [Eqs. (13.23)-(13.25)]. 

The optical effective mass, m*, which, according to Eq. 
(13.26), determines the contribution of the free charge carriers 
to the real part of the dielectric constant, is independent of the 
impurity concentration for a quadratic dispersion relation [Eq. 
(13,29)]. The difference between m* and its value in an intrinsic 
material can be related only to the nonparabolic nature of the 
corresponding band. This effect should be observed in, for ex- 
ample, A!!BY compounds, but obviously not in germanium (or 
Silicon). 


CHAPTER 14 


Heavily Doped Semiconductor with 
an Incompletely Homopolar 


Lattice 


Up to now we have completely neglected lattice oscillations 
in dealing with the problem of the energy spectrum of a semi- 
conductor. Such an approach is reasonable in the case of homo- 
polar crystals like germanium, but when ionic bonds are present 
it raises some doubts. It is well known, for example, that in 
purely ionic crystals consideration ofthe interaction of electrons 
with polarization oscillations of the lattice leads to the concept 
of the polaron [99], a rather radical alteration of the standard 
band model [99-101]. Heavily doped type A''BY. compounds attract 
special interest nowadays. Not being purely ionic, they are, at 
the same time, not completely homopolar; consequently some 
‘‘polaron’’ effects should be expected in them, The question of 
the energy spectrum of a heavily doped semiconductor with the 
interaction of electrons with polarization phonons taken into ac- 
count has been considered in [90]. There the interaction of 
electrons with impurity atoms was completely ignored so that 
the properties of a ‘‘degenerate Fermi gas plus polarization 
phonons’? system were investigated. 

As is known [12], the Green’s function fora system of parti- 
cles which interact with a quantum field can also be written in 
the form of Eq. (5.12),* if by Mis understood the appropriate 


*In [90] causal, not retarded, Green’s functions are used; this, of course, does 
not alter the final results 
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mass operator, One can assume that the contributions of the 
ionic bond in the compounds under consideration is relatively 
small. This makes it possible to use standard perturbation 
theory for calculating the mass operator, In the first approxi~ 
mation the following equation is derived (here and below we are 
following [90]): 


| p—V2mp i: In hv —W, + @ . pe — hv + wo} 


| p+ V2mp | hv — p + Wo Wi AV 4p 


+ Bip, Vv). 


Here A is a constant, Bis a function that does not have a sin-~ 
gularity on the Fermi surface, and  -is the limiting angular 
frequency of the polarization phonons (in energy units). In de- 
riving Eq. (14.1) the phonon energy was assumed to be independ- 
ent of the wave vector and the electron dispersion relation was 
taken in its simplest form (7.13), | 

Near the Fermi surface (p~V2mn) the terms written ex- 
plicitly obviously play a major role in Eq. (14.1) and the function 
B can be neglected. This also justifies the use of the concept of 
‘‘Fermi surface:’? near it the mass operatorisreal and one can 
speak of a dispersion relation. 

The density of states in the case being considered is given 
by Eq. (5.17), Substituting Eq. (14.1) into it, we see that near 
the Fermi surface the function o(£) goes to zero: 


M(p, v) = a In 
(14.1) 


p(E) = ee (14,2) 
1!—aln (j/ =-1) 


Here po is the density of states for an ideal gas, Eq. (5.17); the 
following expression is obtained in [90] for the constant a: 


_ & (0 \* 14.3 
. =a a ( 


where a is the familiar constant for the interaction of electrons 
with polarization phonons, 


Exo and eo being the optical and static dielectric constants, re- 
spectively. 

A result analogous to Eq. (14.2) would also be obtained due 
to the Coulomb interaction of electrons in the first approximation 
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(7.7), Taking screening into account, however, ‘‘smears out?’ 
the logarithmic singularity of the mass operator according to 
Eq. (7.14). In the present case the situation is somewhat differ- 
ent because the field of the polarization phonons is of high 
frequency. On this basis the authors of [90] conclude that 
Screening is absent in the first nonvanishing approximation of 
perturbation theory. It is precisely this circumstance that 
justifies the statement of the problem assumed in [90]: near the 
Fermi surface the interaction of the electrons with polarization 
phonons drastically alters the density of states whereas the 
effect of the impurity in this region, as we have Seen, is in- 
Significant, 

A decrease in the density of states near the Fermi surface 
can apparently explain the phenomenon observed experimentally 
in [102]: the differential conductivity of tunnel diodes fabri- 
cated from group A™!BY compounds has a minimum at low 
voltages (the ascending branch of the V-I characteristic), The 
voltage at which the minimum is observed is proportional to a, 
In this voltage region the tunneling transitions involve electrons 
with energies close to np. 

Equation (14.3) leads to some special features inthe absorp- 
tion of light. Let us consider, for example, a p-type specimen 
and 7<u,. In the assumed approximation electrons in the con- 
duction band behave as if they were almost free. Calculating 
the absorption coefficient from Eq. (13.17) taking into account 
the reality of the mass operator (14.1) one can derive the ex- 
pression 


My 
0] rv on, (“4 ] 
| Vav—A- m, 


ee fe ee 5 
. { a (Wo i. : mM, ( m, hv — A :) 
—_ — | — —— In _ 
27 | By Ms, —-My, Me Mp By, 


Here x must be understood to mean the coupling constant of holes 
with phonons. 

It is seen that because of the polaron effect the absorption 
coefficient does not increase discontinuously from zero for 


m, —-m 
hoo 
hi Vv — \ ~- Ihp an 


m 


but changes relatively smoothly, proportional to x (w,/,)'*, over 
the frequency range. Something similar to this was observed in 
[103]. 
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SUMMARY OF THE RESULTS 
OF CHAPTER 14 


Consideration of the interaction of charge carriers with the 
polarization oscillations of the lattice leads to special charac- 
teristics of a heavily doped semiconductor whose lattice has 
some degree of ionic binding. These special properties emerge 
in the presence of degeneracy, i.e., in the absence of heavy 
compensation. They lead to a sharp decrease in the density of 
states near the Fermi level [Eq. (14.2)]. This leads to a collapse 
of the differential conductivity of tunnel diodes for small forward 
voltages and also to the appearance of a “‘tail’’ in the frequency 
dependence of the band-to-band light absorption coefficient [Eq. 
(14.4)]. 


CHAPTER 15 


Contact of Two Heavily 


Doped Semiconductors 


In a generally macroscopically homogeneous material (all 
we have considered up to now) a number of kinetic effects are 
of interest in an investigation of the energy spectrum. For some 
problems, however, the study of an inhomogeneous system—the 
contact of two heavily doped semiconductors (with the same or 
different types of conductivity)—is essential. Here the energy 
spectra of the two specimens are not as important as are the 
characteristics of the contact itself—first of all, the coordinate 
dependence of the field voltage and chargedensity at the contact. 

In the conventional theory of n-p and n-n’ junctions the cor- 
responding problems are solved by purely classical means. 
The concept of a depletion layer near the contact plays a prom- 
inent role in this case [91]. In the case of heavily doped semi- 
conductors such an approach becomes inadequate for two 
reasons. First, the thickness w of the p-n junction is approxi- 
mately 10-°cm in present-day tunnel diodes; it can be compar- 
able with the characteristic wavelength of the charge carriers, 
which necessitates taking account of quantum corrections. Sec- 
ond, even in the classical case the approximation of a totally 
depleted layer can turn out tobe toosimple because the carriers 
are multiply charged.* 

A direct quantum mechanical calculation of the electric field 
and density of free charge carriers, taking account of their 
interaction with each other, is an exceedingly difficult problem. 


*This circumstance was noted in [104]. 
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For this reason it is advantageous to formulate some variational 
principle for the quantities ofinteresttous, namely, the electric 
field potential U(x) and the Green’s function G(x, x’, f), from 
which the charge density is directly derived.* 

A similar type of variational principle has been formulated 
in [35] and [105] for homogeneous systems, Its extension to the 
spatially inhomogeneous case is given in [106, 107], which we 
will follow in subsequent paragraphs. 

As in [106, 107] we limit ourselves to the Hartree 
approximation, i.e., we take account only of the block diagram 
of Fig. 4 (in an n-p junction the local neutrality condition cer- 
tainly does not exist and the contribution from this graph does 
not become zero by any means). A high field intensity in the 
junction (of the order of V;/w, where V, is the contact potential 
difference) is basic to this. On the other hand, the exchange 
correction to the contact potential difference, according to Eq. 
(10.4), is of the order of A. 

The effects of the ‘‘interaction between bands,’’ as well as 
of the anisotropy, will not play asignificantrole in the problems 
under consideration. Consequently, we can be limited to the 
‘‘one-band’’ approximation of the effective mass operator, as- 
Suming W = p?/2m in this situation. 

Let us first consider a system of charge carriers of one 
Sign (for definiteness electrons). It is convenient to write its 
Hamiltonian in the form 


H =H,+ ef, (15.1) 
where 
a > \ dxa(x, s) (om Vit u) a(x, S), (15.1a) 


and 


gH; == > ( dxdx’ 


a2 940s a(R S)a(x, 3) — 
eo get ce — Xx 


_ £e? \ as = N (x’)a (x, S)a(X, S) - alee ( ax ax n(x) n(x’). 
| \ |x — x’] 


(15.1b) 


€ Jix—x 


S 


*In the present chapter we will utilize the causal Green’s function, G.. The index 
c will be omitted for the most part since the other type of Green’s function is not en- 
countered here. 
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Here n(x) is the localized impurity concentration; g is the 
coupling constant (g =e?/e); s, s’ are spin variables. The term 
with the chemical potential u is added for convenience. 

As is known [12], the shift of the thermodynamic potential 2 
due to the interaction is given by the expression 


g 
AQ = Q(g)—2(0) = f dg’ (Hi (g’)). (15.2) 


Introducing the average value of the operator M,, Eq. (15.1b) 
can be calculated by noting that in the Hartree approximation 


(a(x) a(x’) a(x’)a(x)) = (a(x) a(x) (a(x') a(x’)). (15.3) 
From Eqs. (15.3) and (15.1b) it is easy to obtain 


dx dx’ 


k— x] 


(ff) = Z| [ne (x) — Zn (x)] [Me (x’) — Zn(x’)], (15.4) 


where n,(x) is the local density of electrons with both spin pro- 
jections: 


n, (x) = >; (a(x, s) a(x, S)) = 
=i>) lim G(x, s;x, 5; ¢—¢’) = 2iG(x, x; 0). (19.5) 


5 t’>tto 

The factor of 2 arises here because of the summation over the 
spin: in the absence of spin interactions G(x, s;x’,s’,f—l)= 
b55°G (x, £5 x’, Lf), 

In a locally neutral system the right side of Eq. (15.4) goes 
to zero, just as it should (Chap. 6). 

The equation for the Green function in the approximation 
assumed can be represented in symbolic form: 


AG-' = hGy* + U, (15.6) 
where 
hGy" (x, x',£—-L) = — {in = fet: wh do — x(t —£) 
ot 2m 
(15.7) 
and 
U(x, x1) = U(x) OK — x’) = gO — x A (rex) — Zin’) 


(15.8) 


120 ELECTRON THEORY OF SEMICONDUCTORS 
Obviously, 
yxU (x) = — 4x {n. (x) — Zn(x)} g. (15.8) 


As boundary conditions for Eq. (15.8') we will specifiy the values 
of the potential U on a surface infinitely far away. The Green’s 
function should satisfy the usual type of boundary condition [12]. 

Let us substitute the right side of Eq. (15.4) into Eq. (15.2) 
and integrate by parts, taking Eq. (15.8) into consideration. We 
obtain 


g 


AQ = 8 (Hi) — Vd" (dx U (x) ia nae 
0 


(15.9) 


It is convenient to express the potential energy U(x) in terms of 
the difference of reciprocal Green’s functions with the aid of 
Eq. (15.6), and n,(x) in terms of the Green’s function itself from 
Eq. (15.5). For compactness let us define* the ‘‘spur’’ ofa 
matrix A(x, s; x’, s’,f—t’) by the equation 


SpA = > \ dx lim A(x, s; x’, %; f—’). (15.10) 
5 t’—tto 
We obtain 


g g 
AQ = g (Hl) — ih ( dg'Spg-1 + ih ( dg’SpG> = 
Og’ 6] 
(15.11) 


In light of Eq. (15.7) the last term here is i#SpG, '(G —G,). More- 
over, obviously, SpG—' (aG/ag') = (a/ag’) SpInG. Thus, all inte- 
grals with respect to g’ can be evaluated and Eq, (15.11), taking 
Eqs. (15.4) and (15.5) into account, assumes the form 


AQ = — ihSp {InG —1InG, — G,*(G —G,)} + 
=e ( ss {2iG (x, x, 0) — Zn(x)} {2iG(x’, x’, 0) — Zn(x’)}. (15.12) 


Since the form of the field in the junction will be of greatest 
interest later, it is convenient to introduce the function U(x) ex- 
plicitly in the right side of Eq. (15.12). To this end, we note that 


*This definition differs from the usual one because in Eq. (15.10) there is no in- 
tegration with respect to time. Nevertheless, the left side of Eq. (15.10) possesses 
all the formal characteristics of the usual spur. 
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according to Eqs. (15.4) and (15.8) the expression g<#,> can be 
rewritten in the form 


g <dt> =( dxn. QU) + { dx U(x) Fory v2U (x) — Zn(x)\, 
(15.13) 


Correspondingly, for AQ we obtain, instead of Eq. (15.12), 


AQ = — ihSp {In G —InG, — Gy -(G — G,)} + 
-{- \ ax { ——Uyu — Zn (x) U (x) 4. 2i\ dx U (x)G(x, x; 0) = 
_ ng 
= X[U, G}. (15.12") 


Varying the right side of Eq. (15.12") with respect to U and 
G, we obtain 


Ang = = y2U (x) 4. 4g {2iG(x, x; 0)— Zn(x)}, (15.14) 
x 

Be RG He, Ko SG eS 

2 6G (x’, x; 4 


— U(x’) 8(x — x’). (15.15) 


Obviously the right sides of Eqs. (15.14) and (15.15) go to zero 
if the potential U and the Green’s function G satisfy Eqs. (15.8') 
and (15.7). In other words the expression for AQ is stationary 
with respect to independent variations of U and G around their 
‘‘correct’’ values. 

It is natural to construct the variation AQ from the class of 
functions U’ and G’ which satisfy the same boundary conditions 
as U and G and the related equations of the form (15.6). In other 
words, let us stipulate that G’ is determined from the linear 
Eq. (15.6) with the potential U’. Substituting the solution G’[U’} 
into Eq. (15.12) we obtain a function of U’ only, which should 
also vary. We note that in this case one can also calculate with 
equal success the electron density in terms of the retarded 
Green’s function. It is shown in [106] that under these con- 
ditions the extremum of A®2 is an absolute maximum; in the 
solution of the set of Eqs. (15.14), (15.15) itis unity. 

It should be noted that the variational principle considered 
here has no relationship to the thermodynamic principle of 
minimizing the thermodynamic potential under equilibrium con- 
ditions. It is sufficient to say that for test functions U’ and G’ 
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that are not equal to U and G the functional X[U’, G’] is not equal 
to A®, For this reason the statement concerning the maximizing 
of X should not cause confusion. 

Thus 


DE oe OE ee. RY Gl X (0 G0. (15.16) 
5U - 


in which the equality sign applies only for U’=U and G’=G and 
the functional X is given by Eq. (15.12°). 

The variational principle (15.16) has been formulated expli- 
citly for a system with one type of charge carrier. It is easy to 
extend this to the case of a bipolar system by introducing the 
hole Green’s function from the work of Chaps. 5, 8 [107]. We 
will assume the electron and hole bands are parabolic and neg- 
lect the ‘‘interaction’’ between bands, which apparently plays 
an insignificant role in the problems considered below. Then by 
repeating verbatim all the considerations from above we obtain, 
in the the Hartree approximation, 


AQ = X(U, G., Ga] -= ihSp (Gz) (Ge — Geo) + Gha(Gn — Gro) — 
— (In Ge — In Geo) — (In Gp — In Gao)} — 
— (dxf U (x) VEU (x) —Zatta (x) + Zarta(x)\ — 
P 81g 
— 21 \ axU (x) {Ge (x, x; 0) — G, (x, x; O)}. (15.17) 


Here G, and G, are the Green’s functions for conduction electrons 
and holes; na(x), ‘Na(x), Za and Z, are the localized densities 
and charge (in units of e) of the donors and acceptors, respect- 
ively. The unperturbed reciprocal Green’s functions have the 
form 


ACh te ae gee oe {#5+z - ne ut (x —x’)d(t —1’), 
(15,184) 


RGid (x, x’; {—t)=— {ine + a Vi tha LO — x) dt — Pf), 
(15.18b) 


where m. and m, are the effective masses of the electron and 
hole, and 


Up = —p—A, (15.19) 
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The functional (15.17) is stationary with respect to independ- 
ent variations of G,, Gra and U. This extremum is an absolute 
maximum if the variation is done on the class of functions cited 
above, i.e., in particular, if the test Green functions G, and G, 
are determined from the linear equations 


AG. =hG,' + U’, (15.20a) 
and 
AG, = RG, + U', (15.20b) 


by then varying X with respect to U’. 

The variational principle in the form of (15.16), (15.12) and 
(15.17) can be used directly for an approximate calculation of 
the potential and charge density at the contact or in any other 
highly inhomogeneous system. In this case any quantum or 
classical effects which are simply associated with the effect of 
free charge carriers on the contact characteristics can be 
studied. 

Let us consider first of all the simplified one-dimensional 
problem of the contact of two heavily doped semiconductors with 
one type of conductivity (e.g., a degenerate n—n’ junction; 
T=0). 

We will assume the junction is ‘‘abrupt’’ in the normal 
technical sense of the word: 


—n= 0, 

nj n=const, *< (15.21) 

= n’ = const, x > QO, 

for which n>n’., We choose a test function in the form 
U (x) = V, : exp { (1 26) <0 
=I -'- 0 ) | w |’ (15.22) 
U(x) Vell — exp| — 2(1-+ 8). x>0. 
1-6 wf 


Here 6 and w are parameters to be determined from the vari- 
ational principle. It is easy to see that w is nothing more than 
the junction width, and 6 characterizes its.asymmetry (when 
6 -1, the field intensity F == —dU/edxis symmetrical with re- 
spect to the point x = 0). 

Under the condition that } is small the Green’s function in 
the test field (15.22) can be calculated by semiclassical methods 
(as in Chap. 8). Substituting Eq. (15.22) into Eq. (8.17), calculat- 
ing n-(x) and carrying out the variation, we obtain equations for 
«and 6. Thus we obtain for the junction width 
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w? = wo(l — Q), (15.23) 
where w, is the classical junction width. An explicit expression 


for this can be found in [106]. The quantity Q is a quantum 
correction term: 


Oi ee, (15.24) 
m (4 + €)(1+ 2) (2 +€) 
Here* 
V 3 eee 
aV/i-t- a, oe nar 
UW 1+ 5 


and 6, represents the classical value of 6; the function 6, (V2/) 
is tabulated in [106]. As V,/v varies from 0 to 1 the coeffi- 
cient of A’ in Eq. (15.24) changes monotonically from 1/3 to 7/2. 
It is seen that the quantum correction reduces the junction width 
to some extent. 

The effect of free charges on the charging capacitance C of 
the degenerate p-n junction has been investigated in [107] 
with the aid of the variational principle (15.17). For sufficiently 
small contact potential differences the result can differ by a con- 
siderable factor from that given bynormal depletion layer theory. 
That is, for T<w one obtains 


ees (15.25) 


where C, is the capacitance calculated in the depletion layer 
approximation. 


SUMMARY OF THE RESULTS 
OF CHAPTER 15 


The contact of two degenerate semiconductors (with the same 
or different types of conductivity) has special properties as- 
sociated with quantum effects (in a thin junction) and with the 
effect of free charge carriers. For sufficiently small contact 
potential differences the latter significantly reduces the junction 
capacitance, compared with the usual result of depletion layer 
theory [Eg. (15.25)}. Quantum corrections affect the junction 
width (reducing it somewhat); however, their role has not yet 
been investigated thoroughly. 


*In the case under consideration it is obvious that V, <u. 


to. = 
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